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Abstract 

Over  the  last  decade,  quantum  information  experiments  with  trapped  ions  have  demon¬ 
strated  essential  steps  towards  quantum  computing  and  quantum  simulation.  Large  fields 
are  required  to  achieve  strong  coupling  to  the  ions  via  dipolar  interactions,  and  so  we  fab¬ 
ricated  transmission  line  microresonators-capable  of  producing  large  fields  in  a  standing 
wave  at  resonance— for  eventual  integration  into  2D  ion  trap  structures.  The  resonators 
were  superconducting  to  minimize  loss  and  maximize  quality  factor.  We  fabricated  the 
resonators  as  two  dimensional  coplanar  waveguides  in  niobium  on  R-plane  sapphire  using 
optical  lithography.  Resist  was  patterned  on  the  niobium  using  optical  lithography,  devel¬ 
oped,  then  reactive-ion  etched  to  transfer  the  pattern  into  the  niobium.  The  resonators  were 
cooled  and  tested  in  a  cryogenic  probe  station  and  characterized  with  a  network  analyzer. 
Additionally,  the  resonator  geometry  was  reproduced  in  commercial  microwave  simulation 
software.  Results  from  our  fabricated  resonators  showed  first-resonance  quality  factors  of 
1.2x10“^  at  3.23  GHz  at  device  temperatures  of  3-4  K. 
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Chapter  1 


Introduction 


Resonators  are  one  of  the  most  fundamental  types  of  passive  circuit,  and  though  they  have 
been  around  for  over  a  century,  current  research  is  still  finding  new  uses  for  them.  By 
scaling  the  size  down  using  microfabrication  techniques,  it  is  possible  to  build  resonators 
that  operate  in  the  microwave  regime,  extending  their  applications  beyond  the  realm  of  their 
traditional  place  in  lumped-element  analog  electronics.  Recent  applications  include  use  in 
microwave  kinetic-inductance  detectors [29],  multiplexed  SQUIDs[10],  and  most  pertinent 
to  the  work  discussed  in  this  thesis,  ion  trapping  and  quantum  information  processing. 

1.1  Quantum  Information  Processing  and  Ion  Trapping 

Quantum  information  processing  is  a  relatively  new  field  that  seeks  to  expand  the  realm 
of  traditional  computing  by  interfacing  with  quantum  mechanical  phenomena  to  build  and 
measure  quantum  bits,  or  qubits.  Algorithms  already  exist  within  the  field  that  offer  sig¬ 
nificant  speedup  to  perform  traditionally  difficult  tasks  such  as  factoring  large  primes  or 
searching  large  databases.  A  particularly  promising  way  of  building  qubits  for  quantum 
computation  comes  in  the  field  of  ion  trapping. 

1.1.1  Ion  trapping 

Ion  trapping  comes  in  many  forms,  but  recently  there  has  been  significant  development  in 
the  area  of  planar  ion  trapping.  A  planar  ion  trap  is  functionally  identical  to  a  standard 
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linear  RF  Paul  trap[17],  in  that  it  creates  an  oscillating  quadrature  electric  field  with  a  local 
minima  able  to  capture  ions  in  two  dimensions  [18].  The  difference  lies  in  the  geometry  of 
the  trap,  which  adapts  the  confinement  electrodes  from  their  standard  three-dimensional 
configuration  into  a  two-dimensional  planar  structure. 

Planar  RF  ion  traps  have  many  advantages  over  their  3D  cousins.  Most  importantly, 
their  planar  structure  allows  for  rapid  prototyping  with  standard  microfabrication  tech¬ 
niques.  Additionally,  they  hold  promise  as  a  scalable  architecture,  since  any  one  trap  may 
hold  several  ions,  and  multiple  traps  may  be  fabricated  over  a  small  area  and  potentially 
used  as  a  network  of  processing  and  memory  zones  [23].  Lastly,  they  offer  advantages  in 
cryogenic  ion  trapping,  such  as  the  ion  trapping  currently  being  researched  in  MIT’s  Center 
for  Ultracold  Atoms  (CUA). 

1.1.2  Current  ion  trapping  research  at  the  CUA 

We  have  previously  demonstrated  the  ability  to  trap  ®®Sr+  in  the  CUA  using  a  cryogenic 
planar  ion  trap[2][3][12]  with  lifetimes  of  over  2500  min.  Operation  at  cryogenic  temper¬ 
atures  is  particularly  useful  in  ion  trapping  because  of  the  high  vacuum  it  generates-ion 
lifetime  is  typically  cut  short  through  interactions  with  residual  gasses.  The  trap  was  planar 
both  to  simplify  fabrication,  and  to  allow  contact  of  the  trap  substrate  (and  subsequently, 
electrodes)  directly  against  the  cooling  stage  of  the  cryostat  in  which  the  apparatus  is  tested. 
The  strontium  ions  were  laser  ablated  from  a  SrTiOs  target,  then  Doppler  and  sideband 
cooled  into  the  motional  ground  state  for  trapping.  Once  trapped,  the  ions  were  detected 
with  an  optical  setup  that  observed  their  fluorescence. 

1.2  Towards  an  Ion  Trap  with  an  Integrated  Superconduct¬ 
ing  Resonator 

Though  trapping  has  already  been  demonstrated,  we  now  seek  to  interact  with  the  trapped 
particles  and  thus  take  a  step  towards  quantum  processing.  Specifically,  we  hope  to  trap 
SrCl2  molecules  and  excite  their  long-lived  rotational  states[l].  These  molecules  have  a 
transition  near  6  GHz,  and  so  we  should  be  able  to  demonstrate  coupling  by  placing  them 
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in  a  strong  electric  field  oscillating  at  that  transition  frequency. 

Rather  than  having  to  purchase  a  military-grade  microwave  source  capable  of  generating 
hundreds  of  watts  of  microwave  power  to  achieve  this  coupling,  a  transmission  line  cavity 
resonator  can  be  used  to  accumulate  microwave  photons.  Operated  at  resonance,  the  pho¬ 
tons  in  the  cavity  form  a  standing  wave  with  a  high  electric  field  (orders  of  magnitude  larger 
than  the  source  wave)  at  the  field  antinodes.  By  integrating  a  transmission  line  resonator 
into  the  planar  trap  geometry,  we  hope  to  be  able  to  maneuver  the  ion  near  the  localized 
nodes  of  high  electric  field  and  interact  with  it  in  the  strong  coupling  regime. 

1.2.1  Designing  the  resonator  for  our  needs 

As  a  building  block  towards  the  larger  goal  of  quantum  computation,  this  thesis  is  focused 
on  the  fabrication  of  suitable  resonators  for  trap  integration.  Precisely  how  many  photons 
are  able  to  be  trapped  in  the  resonator  cavity-and  thus  how  high  an  electric  field  can  be 
achieved-is  determined  by  a  critical  parameter  called  the  quality  factor,  or  Q.  We  desired 
a  high  quality  factor,  and  so  the  transmission  line  resonators  described  in  this  thesis  were 
designed  using  a  superconducting  coplanar  waveguide  to  minimize  energy  lost  through  dis¬ 
sipation  in  the  metal  and  dielectric.  The  effects  of  superconductivity  on  the  resonator  as  a 
whole  and  the  quality  factor  are  the  subject  of  discussion  in  Chapter  3. 

In  Chapter  2  we  describe  the  complete  fabrication  of  superconducting  niobium  resonators 
designed  to  achieve  very  high  quality  factors  with  a  first  resonance  at  3  GHz.  The  resonators 
are  completely  planar,  and  the  process  involves  only  the  manipulation  of  a  single  niobium 
layer  on  a  sapphire  substrate.  Briefly,  the  conductor  layer  is  deposited,  photoresist  is 
patterned  on  its  surface,  and  then  the  resist  pattern  is  transferred  into  the  conductor  by 
means  of  reactive-ion  etching. 

The  majority  of  the  resonators  fabricated  were  characterized  using  a  cryogenic  probe 
station  and  network  analyzer.  Information  on  the  probe  station,  other  apparatus,  and 
procedures  may  be  found  in  Chapter  5.  Lastly,  the  results  of  of  our  work  are  explored  in 
Chapter  6. 
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Chapter  2 


Materials  and  Fabrication 


This  chapter  details  the  processes  developed  in  the  production  of  our  superconducting  res¬ 
onators.  Much  of  the  previous  work  done  examining  these  devices  was  done  with  either 
aluminum  or  niobium  as  the  superconductor,  and  silicon  or  sapphire  as  the  dielectric  sub¬ 
strate.  Many  factors  affect  the  higher-order  properties  of  superconducting  resonators  (such 
as  kinetic  inductance),  but  in  terms  of  the  first-order  lumped-element  model,  the  resonators 
may  be  fabricated  out  of  any  conductor  and  substrate  combination.  Since  we  began  this 
investigation  with  no  pre-existing  devices  to  adapt,  we  opted  to  use  materials  which  have 
both  already  been  shown  to  produce  high  quality  factors,  and  which  we  were  able  to  fab¬ 
ricate  from  day  one:  niobium  and  sapphire.  Since  the  expertise  for  depositing  niobium 
films  already  existed  within  the  group,  we  used  the  existing  niobium  process  and  built  the 
remainder  of  our  fabrication  steps  on  top  of  that. 

2.1  The  Material  Stack 

Fabrication  began  with  a  330pm  thick,  5  cm  crystalline  R-plane  sapphire  wafer.  The  wafers 
were  ordered  from  Kyocera  Industrial  Ceramics,  and  they  were  epitaxially  polished  on  one 
side,  and  chemically  polished  on  the  other.  Since  the  wafers  were  R-plane,  they  exhibited 
dielectric  anisotropy  between  the  two  axes  parallel  with  the  cut  of  the  wafer.  These  two  axes 
were  the  A-axis  and  the  C-axis,  which  according  to  the  Kyocera  data  sheet  had  effective 
dielectric  constants  of  9.3  and  11.5  respectively.  For  our  devices,  all  of  our  calculations  were 
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done  with  regard  to  the  C-axis.  The  reason  for  choosing  the  higher  dielectric  constant  was 
to  force  more  of  the  resonating  field  out  of  the  sapphire,  where  the  high  dielectric  constant 
made  field  density  energetically  expensive,  and  put  more  of  the  field  into  the  surrounding 
vacuum,  where  it  could  interact  with  the  ion. 

Niobium  was  used  primarily  because  the  deposition  capability  for  this  material  already 
existed  within  the  lab.  We  additionally  tried  using  niobium  nitride,  but  NbN  is  highly 
variable  in  material  phase  and  Tc[13][16],  based  on  the  temperature  and  gas  ratio  with 
which  it  is  sputtered.  Despite  the  higher  critical  temperature  of  NbN,  we  decided  that  the 
simpler  material  of  niobium  was  the  best  route  in  the  early  stages  of  developing  a  high 
quality  factor.  Previous  literature  suggested  niobium  was  a  suitable  superconductor  for 
quality  factors  in  the  range  of  10®  -  lO'^  [9]  [11]  [5]  [22],  and  our  probe  station  tests  operated 
at  around  3-4  K,  well  below  the  critical  temperature  of  niobium.  Additionally,  we  were 
advised  by  Dr.  Ben  Mazin  not  to  bother  with  NbN,  and  that  other  issues  such  as  sloping 
sidewalls  were  more  critical  to  the  resonator  fabrication. 

All  of  the  primary  deposition  and  fabrication  was  realized  on  the  epi-polished  side  of 
the  wafer,  although  the  backside  was  also  coated  for  purposes  of  heat  conduction  and  to 
provide  a  ground  layer.  Niobiumwas  deposited  on  the  wafer  by  means  of  an  AJA  sputtering 
system^  located  in  the  Quantum  Nanostructures  and  Nanofabrication  Laboratory  at  MIT 
and  shown  in  Figure  2-1.  In  the  following  subsections  we  describe  the  procedure  used  to 
deposit  niobium  in  the  AJA  sputtering  system  and  build  the  Nb-sapphire-Nb  “material 
stack.” 

2.2  Resonator  Fabrication 

The  process  began  by  first  coating  the  backside  of  the  wafer  with  a  thin  film  of  sputtered 
niobium  from  a  99.95%  pure  target.  This  layer  was  never  patterned,  and  was  only  necessary 
to  help  the  wafer  absorb  thermal  radiation  from  the  hot  chuck  when  the  polished  side  was 
being  coated.  In  the  following  paragraphs,  the  backside  coating  procedure  is  neglected  for 
the  more  complete  frontside  coating  procedure,  but  is  detailed  afterwards.  A  wafer  cleaning 
^Model  ATC  2000  thin  film  deposition  system,  built  by  AJA  International,  Inc 
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Figure  2-1;  The  QNN  Laboratory’s  AJA  Sputterer.  (A)  denotes  the  sputtering  system’s 
load- lock  entry  point.  (B)  denotes  the  sputtering  chamber,  30  cm  in  diameter  and  40  cm 
in  height. 
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procedure  was  omitted  because  the  front  surface  arrives  from  the  factory  epitaxially  smooth. 
We  rely  on  the  heat  and  low  pressure  of  the  sputtering  deposition  chamber  to  remove  any 
water  contamination  deposited  on  the  wafer  during  transport  from  the  cleanroom  to  the 
AJA. 

The  wafer  was  loaded,  frontside-down,  into  an  inconel  chuck  cover  and  was  screwed 
into  the  base  of  the  chuck.  The  chuck  and  wafer,  now  stuck  together,  were  placed  into  the 
load  lock  of  the  sputtering  system.  The  load-lock  turbo  was  engaged  and  once  its  pressure 
reaches  3xl0~^  Torr,  the  main  chamber  gate  wasopened  and  a  magnetic  lever  arm  was  used 
to  push  the  sample  into  the  main  chamber  and  load  it  onto  the  main  chamber  chuck-mount. 
Once  the  transfer  arm  was  out  of  the  chamber  and  the  gate  was  closed,  the  chuck  holding 
the  sample  was  set  to  rotate  at  50  rpm  and  the  relative  height  level  was  adjusted  to  47  mm 
for  the  remainder  of  the  sputtering  procedure.  This  height  setting  of  47  mm  corresponded 
to  a  target-sample  distance  of  18  cm. 

The  temperature  control  was  set  to  850°  C,  and  after  7-8  min,  the  chuck  and  sample 
reached  temperature.  We  waited  an  additional  5  min  after  first  reaching  850°  C  to  stabilize 
the  temperature  and  also  allow  any  outgassing  contaminants  to  dissipate.  During  this 
time  period  the  pressure  typically  remained  below  1  x  10“®  Torr.  Pressures  above  this  level 
indicated  a  leak  in  the  system  or  large  amounts  of  outgassing  contaminants.  Ideally,  when 
the  wafer  had  finished  outgassing  the  pressure  baseline  was  in  the  low  1x10“®  Torr  range. 

Now  we  were  ready  to  perform  the  pre-sputtering  maintenance.  Keeping  the  argon 
cylinder  regulator  closed,  we  set  the  argon  flow  rate  into  the  main  chamber  to  100  seem 
and  waited  until  the  flow  rate  dropped  to  below  1  seem,  in  order  to  drain  any  argon  or 
contaminants  trapped  in  the  flow  line.  We  then  flushed  the  gas  lines  by  opening  the  argon 
cylinder  regulator,  restoring  the  100  seem  flow. 

After  flushing  the  lines,  we  performed  a  short  target  cleaning  procedure  to  sputter  away 
the  top  layer  of  the  niobium  target,  which  may  have  been  contaminated  from  other  users, 
venting,  or  slow  leaks  in  the  chamber.  While  cleaning  the  target,  the  shutter  on  the  target 
remained  closed,  so  that  any  particulates  ejected  from  the  contaminated  surface  of  the  target 
would  be  collected  on  the  inside  of  the  shutter  and  not  spoil  the  surface  of  the  sample.  The 
power  for  the  sputtering  target  was  sourced  from  a  DC  gun  at  a  maximum  possible  current 
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of  1  A.  With  the  argon  flowing  at  100  seem,  the  DC  gun  current  was  set  to  30%  (300  mA), 
and  the  chamber  pressure  was  set  at  30  mTorr.  The  pressure  was  immediately  lowered  to 
3  mTorr  once  sparked,  and  the  target  plasma  was  left  on  to  clean  for  3  min. 

Using  the  parameters  of  the  last  paragraph  in  the  sputtering  system,  niobium  accumu¬ 
lated  on  the  target  at  27A/min,  measured  by  system’s  built-in  quartz  oscillator  thickness 
monitor.  To  achieve  the  thick  (320  nm)  niobium  we  used,  the  target  shutter  stayed  open 
for  2  hours,  after  which  the  target  shutter  was  closed,  the  plasma  was  shut  off,  and  the  gas 
flow  was  set  to  zero.  For  sample  cooldown  the  cryo  pump  was  fully  opened  to  bring  the 
pressure  down  as  much  as  possible,  and  the  sample  was  cooled  for  at  least  two  hours,  al¬ 
though  typically  it  remained  in  the  chamber  overnight.  After  sputtering,  the  substrate  was 
immediately  placed  into  a  nitrogen  drybox  until  the  next  step  in  processing,  to  minimize 
oxidation. 

The  backside  sputtering  procedure  was  identical  to  the  frontside,  except  that  the  heating 
was  omitted  and  the  sputtering  time  was  reduced  from  2  hours  to  5  min,  just  enough  to  give 
a  visible  layer.  Note  that  the  backside  sputtering  takes  place  first,  since  it  was  important 
for  the  even  heating  of  the  wafer  during  frontside  deposition. 

2.2.1  Photolithography  and  etching 

Once  the  two  niobium  layers  were  deposited,  the  wafer  was  ready  for  processing.  To  fabricate 
our  resonators,  micron-level  precision  was  needed,  and  so  processing  was  done  using  an 
ultraviolet  optical  lithography  setup  and  reactive  ion  etching.  A  schematic  showing  the 
layout  of  the  resonator  was  shown  in  Figure  2-2.  The  wafer  mask  was  designed  in  the  GDSII 
format  using  the  program  LayoutEditor^,  and  manufactured  at  Advanced  Reproductions 
Corp.^  on  a  10  cm  square  soda- lime  mask  patterned  with  1.5  mm  chrome.  Because  the 
photoresist  used  was  negative,  the  mask  was  a  negative  mask,  which  means  that  all  the 
electrodes  of  the  fabricated  devices  were  blank  areas  on  the  mask.  Correspondingly,  gaps 
on  the  final  resonator  had  their  chromium  equivalents  on  the  mask.  Once  the  mask  was 
designed  and  fabricated,  the  remainder  of  the  processing  was  performed  in  the  NSL  clean 

^http;/ /www.layouteditor.net/ 

^http://www.ad  vancerepro.com 
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room.  All  of  the  following  processes  except  baking  occured  at  room  temperature,  which  for 


the  NSL  clean  room  was  23±1°  C. 


Figure  2-2:  Schematic  of  the  resonator  showing  the  contact  pads  (grey  electrodes)  and 
resonating  center  electrode  (red).  Not  to  scale. 


We  began  by  spinning  NR9-3000P  photoresist^  onto  the  wafer  at  3000  rpm  for  40  sec, 

yielding  a  resist  thickness  of  approximately  three  micrometers  as  seen  in  SEM  images.  Next, 

the  wafer  was  prebaked  on  a  hot  plate  at  90°  C  for  5  min  after  which  it  was  moved  to  the 

Tamarack  exposure  system  The  Tamarack  lamp  was  turned  on  and  allowed  to  warm  up 

for  10  min,  after  which  the  exposure  intensitywas  measured.  The  typical  exposure  intensity 

was  3200  /iW/cm^,  and  at  that  intensity  the  resist  was  exposed  for  90  sec,  yielding  a  total 

exposure  of  288  mJ/cm^.  The  NR9-3000P  SOP  suggests  much  hotter  bakes  (150°  C)  than 

we  were  using,  but  we  devised  this  recipe  to  avoid  heating  the  niobium  too  much  and 

possibly  increasing  the  thickness  of  a  loss-inducing  oxide  surface  layer.  Additionally,  the 

SOP  states  that  only  an  exposure  of  21  mJ/cm^  was  necessary,  but  early  tests  revealed  the 

change  in  bake  temperatures  also  required  much  longer  exposures. 

‘‘From  Futurrex  http://www.futurrex.com/ 

“Model  PRX-200-6  (Hg  Lamp),  Tamarack  Scientific  Co. 
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The  mask  was  manually  placed  onto  the  sample,  and  a  vacuum  seal  was  placed  around 
the  mask  to  pump  out  any  air  between  the  resist  and  the  mask  and  ensure  intimate  contact. 
Typically  the  interference  fringes  seen  in  the  yellow  clean  room  light  were  limited  to  two 
or  three  localized  bands,  indicating  some  small  grains  of  dust  were  present  on  the  surface 
but  that  intimate  contact  to  within  1-2  pm  was  achieved  over  most  of  the  surface.  The 
mask-wafer  pair  was  then  moved  directly  under  the  Tamarack  lamp  and  was  exposed  foi 
90  sec.  Following  the  exposure,  the  wafer  was  baked  at  90°  C,  this  time  for  2  min. 

Now  the  sample  was  ready  for  development  and  etching.  The  sample  was  submerged  in 
room-temperature  RD6®  and  agitated  lightly  in  the  developer  for  16  sec  after  which  it  was 
immediately  rinsed  by  placing  it  under  the  flow  of  DI  water  for  30  sec.  After  rinsing,  the 
sample  was  blown  dry  with  nitrogen  and  taken  to  the  optical  microscope  for  examination^, 
the  results  of  which  can  be  seen  in  Figure  2-3. 


100um 


Figure  2-3:  Sample  amccO  10.04  4-finger  coupling  capacitor,  patterned  into  NR9-3000P  resist 


If  features  were  not  fully  developed,  the  sample  was  taken  back  to  the  chemical  bench 
and  developed  for  3  more  seconds,  following  the  same  development  procedure  as  before.  The 
criteria  for  “fully  developed”  was  only  that  the  smallest  feature  gaps  (in  our  case,  10pm) 
were  visible  in  the  microscope.  Once  features  were  satisfactorily  resolved  in  the  resist,  the 
sample  was  baked  one  last  time  at  90°  C  for  2  min  and  could  then  be  etched.  Etching  was 

®Froni  Futurrex  http://www.futurrex.com/ 

’’It  should  noted  that  the  microscope  must  be  fitted  with  a  UV  filter  before  examination,  lest  it  further 
expose  the  portion  of  the  sample  being  viewed. 
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performed  in  a  Plasma  Therm  790  series  RIE.  The  process  begins  by  cleaning  the  chamber 
for  5  min  using  the  preconfigured  “NEWCLEAN”  recipe,  which  scrubs  the  inside  of  the 
chamber  with  the  same  set  of  gases  that  the  etching  process  used.  Afterwards,  the  sample 
was  placed  onto  the  conductive  center  of  the  RIE  chamber.  The  lid  was  closed,  and  once 
the  vacuum  seal  was  made  the  CF4/O2  process  began.  The  chamber  was  reduced  to  a 
base  pressure  of  8x10“^  Torr  and  held  at  that  level  for  10  sec.  Next,  to  flush  the  lines 
and  chamber,  a  gas  mixture  flowed  for  1  min  at  20  mTorr,  with  a  ratio  of  1.5  seem  O2 
to  15  seem  CF4.  After  the  lines  were  flushed,  the  plasma  was  sparked  and  held  at  150  W 
of  RF  power  and  the  sample  was  etched  for  24  min.  Finally,  the  chamber  was  pumped 
down  to  2x10“^  Torr,  held  there  for  10  sec,  and  vented  to  atmosphere.  A  low-zoom  visual 
inspection  of  the  surface  took  place  using  the  Dektak’s  camera  to  ensure  that  the  niobium 
had  been  fully  etched  away  in  the  exposed  areas.  The  Dektak  camera  was  used  since  it 
zooms  in  on  the  sample  from  a  45°angle  while  the  sample  was  illuminated  from  above.  This 
made  it  easier  to  spot  any  thin  residual  layers  of  niobium  that  were  not  fully  etched  away. 

If  niobium  remained,  a  second,  shorter  round  of  etching  in  the  RIE  took  place,  using 
the  exact  same  process  outlined  above  except  with  1-2  min  of  etching  time.  Otherwise, 
if  the  etching  was  satisfactory,  the  sample  was  cleaned  of  the  remaining  photoresist  by 
sonicating  it  in  an  acetone  bath  for  5  min.  Normally  the  photoresist  could  be  removed 
just  by  submerging  it  in  acetone,  but  the  etching  procedure  interacts  with  the  resist  and 
the  resulting  photoresist  post-etch  was  much  more  difficult  to  remove,  hence  the  sonication. 
Lastly,  the  sample  was  gently  rinsed  with  acetone,  methanol,  IPA,  and  then  deionized  water 
in  that  order,  after  which  it  was  blown  dry  with  nitrogen.  The  final  result  was  visible  in 
Figure  2-4. 

2.2.2  Gold  Contact  Pads  and  liftoff 

Although  the  resonator  was  fundamentally  complete  after  etching,  in  order  to  test  it  in  our 
probe  station  we  needed  to  add  a  gold  layer  to  the  surface  of  the  contact  pads,  so  that 
the  probes  could  make  good  electrical  contact.  Initially,  we  probed  the  niobium  while  it 
was  bare.  This  approach  did  not  work  well,  because  in  the  cryogenic  regime  the  niobium 
was  very  hard,  and  the  probes  had  a  difficult  time  making  contact.  Even  small  vibrations 
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Figure  2-4:  Optical  microscope  images  of  the  etched  resonator  from  sample  amccOOl-01 

would  move  and  disturb  the  contacted  probes  enough  to  ruin  transmission  measurements. 
Additionally  the  unreliable  electrical  contact  made  it  difficult  get  consistent  results-without 
gold  contact  pads,  probing  and  reprobing  the  same  resonator  yielded  wildly  varying  (±10  dB 
across  the  entire  spectrum)  transmission  characteristics. 

The  gold  contact  pads  were  added  to  the  device  using  a  liftoff  procedure  very  similar 
to  the  etching  procedure,  since  the  liftoff  resist  used-NR9-3000PY-was  a  variant  of  the 
photoresist  used  in  etching,  NR9-3000P.  First,  the  gold  contact  mask  was  made  by  editing 
the  device  mask  in  LayoutEditor  and  sending  it  to  be  printed  on  an  emulsion  transparency* 
.  The  transparency  had  very  rough  features,  with  a  minimum  feature  size  of  25pm,  but 
since  the  devicewas  already  accurately  patterned,  the  gold  pattern  on  top  could  be  very 
rough. 

The  resist  was  spun  onto  the  patterned  device  at  3000  rpm  for  40  sec  and  then  prebaked 
at  90°  C  for  2  min.  The  samples  were  then  wrapped  in  foil,  taken  out  of  the  NSL,  and 
brought  to  the  MJ83  broadband  UV  mask  aligner  in  the  Microsystems  Technology  Labora¬ 
tory  (MTL).  At  the  mask  aligner,  the  transparency  mask  was  taped  to  a  10  cm  glass  plate 
and  the  pattern  was  aligned  with  the  resist-coated  wafer.  Once  aligned,  hard  contact  was 
forced  between  the  mask  and  wafer,  and  the  combination  was  exposed  with  the  internal 
mercury  lamp  for  2  min.  Just  as  in  the  etching  process,  the  exposed  wafer  was  then  baked 
at  90°  C  for  2  min  and  developed  in  RD6  for  16  sec.  In  this  case,  however,  the  resistwas 
*5080  dpi  transparency  printed  from  Pageworks,  Inc. 
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not  post-baked  after  development,  since  SEM  examination  showed  that  the  NR9-3000PY 
reflowed  if  post-baked,  changing  its  sidewall  profile  from  an  undercut  to  a  shallow  slope  in¬ 
appropriate  for  liftoff.  After  development,  the  wafer  was  mounted  onto  a  fixture  and  placed 
into  an  evaporator,  where  the  wafer  collected  10  nm  of  titanium  (for  adhesion)  followed  by 
60  nm  of  gold.  Once  the  gold  was  deposited,  the  wafer  was  then  sonicated  in  acetone  for 
5  min  to  remove  the  remaining  resist  and  with  it,  the  excess  gold.  The  sonication  was  not 
typically  necessary  for  this  type  of  liftoff,  but  was  performed  to  make  absolutely  sure  no 
hardened  resist  was  leftover  from  the  etching  procedure.  It  was  finally  cleaned  in  acetone, 
followed  by  methanol  and  isopropanol,  and  then  rinsed  in  dionized  water  and  blown  dry 
with  nitrogen.  The  final  result  can  be  seen  in  Figure  2-5. 


Figure  2-5:  The  fully  fabricated  set  of  straightline  resonators  mounted  in  a  copper  frame 


Chapter  3 


Theory  of  Superconducting 
Resonators 


To  accurately  understand  the  measurements  we  made,  we  examined  the  theory  behind 
superconducting  resonators.  Starting  with  the  basic  models  of  transmission  line  theory  and 
lumped  element  resonators,  we  were  able  to  interpret  our  results  and  gain  an  intuition  for 
how  to  improve  the  devices.  The  following  sections  set  up  the  theoretical  understanding 
necessary  to  analyze  the  results  presented  in  Chapter  6. 

3.1  Transmission  Line  Model 

Everything  about  the  resonator  can  be  analyzed  using  standard  transmission  line  theory. 
In  this  section,  we  break  down  the  resonator  first  as  a  transmission  line,  then  further  into 
the  constituent  parameters  of  the  transmission  line’s  distributed  circuit  model  in  order  to 
understand  how  the  physical  shape  of  the  device  affects  its  operation. 

3.1.1  The  Transmission  Line  as  a  Circuit  Model 

Though  the  parallel  RLC  circuit  gives  us  an  excellent  description  of  the  overall  circuit 
parameters  and  intuition  as  to  what  kind  of  behavior  to  expect,  it  paints  an  incomplete 
picture.  In  a  transmission  line,  the  phase  and  magnitude  along  the  length  of  the  resonator 
can  be  most  fundamentally  examined  by  breaking  the  system  down  into  a  distributed  cir- 
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Figure  3-1:  Segment  /S.z  of  a  general  two-conductor  transmission  line  represented  in  its 
distributed  circuit  model  form. 


cuit  model,  where  an  infinitesimal  length  of  the  TL  is  represented  by  the  circuit  model  of 
Figure  3-1.  This  formulation  is  more  thoroughly  explained  in  the  literature  [21].  We  will 
simply  note  that  a  uniform  lossless  transmission  line  with  characteristic  impedance  Zq  with 
propagation  constant  /3  has  the  standard  solution 


F(z)  =  Vb[e"^^^  +  re^'^^] 


I{z)  =  -  Fe^^"]. 

Zo 


where  F  is  the  voltage  reflection  coefficient,  and  V{z,t)  —  Re{V {z)e^‘'^^} .  The  device  under 
consideration  is  most  accurately  described  by  an  open-circuit  half-lambda  (A/2)  resonator, 
for  which  F  =  1.  The  input  impedance  of  such  a  lossless  resonator  is 


Zin^  Zocoth{jl3)l  (3.3) 

The  resonant  modes  occur  whenever  the  length  I  is  a  multiple  of  A/2,  where  A  =  2%/^. 
The  resulting  voltage  distribution  V{z)  for  the  first  3  resonant  modes  is  seen  in  Figure  3-2. 

3.1.2  Physical  Geometry  of  the  Device 

Since  the  transmission  line  model  is  generalized,  the  device  can  be  designed  in  any  particular 
geometry  and  still  maintain  the  fundamental  properties  of  the  TL  resonator.  As  long  as  the 
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Figure  3-3;  Side  view  of  the  CPW  geometry  showing  the  characteristic  parameters  width 
of  center  signal  line  w,  spacing  between  signal  line  and  ground  s,  and  thickness  of  the 
conductors  t. 


Figure  3-4;  Top  view  of  the  fabricated  CPW  line  showing  the  resonator  transmission  line 
(left  electrode)  and  input  transmission  line.  The  gap  is  the  coupling  capacitor. 
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parameters  R,  L,  G,  and  C  per  unit  length  can  be  found,  the  general  model  is  useful.  The 
device  addressed  in  this  report  has  been  designed  as  a  grounded  coplanar  waveguide  (CPW) 
transmission  line  resonator,  as  shown  from  above  in  Figure  3-4  and  whose  cross-section  can 
be  seen  in  Figure  3-3.  The  grounded  CPW  consists  of  a  thin  flat  dielectric  sandwiched 
on  either  side  by  a  conductor,  one  side  of  which  is  etched  to  form  a  signal  line  bounded 
on  either  side  by  ground  planes.  In  our  particular  case,  we  chose  to  use  niobium  for  the 
conductor  and  sapphire  for  the  dielectric.  The  niobium  was  chosen  for  its  superconducting 
properties  and  convenience  of  fabrication.  The  sapphire  makes  an  excellent  dielectric  due 
to  its  low  loss  tangent  and  high  permittivity  as  discussed  in  Chapter  2. 

With  the  Nb-sapphire  material  stack,  the  conductance  G  is  negligible  and  can  be  re¬ 
moved  from  the  transmission  line  circuit  model.  However,  there  is  no  simple  field  solution 
to  solve  like  that  which  exists  for  a  parallel  plate  transmission  line.  Fortunately,  previous 
work  [24]  [8]  has  been  done  to  find  a  conformal  mapping  technique  for  the  CPW  that  yields 
the  geometric  contribution  to  inductance  and  capacitance  per  unit  length 


rm  _  /^O 

‘  ”  4  K{ko) 


(3.4) 


In  these  equations  K  is  the  complete  elliptic  integral  of  the  first  kind  with  the  geometric 
arguments 


ko  = 


w 

w  +  2s 


(3.6) 


A:'  =  (3.7) 

Though  these  equations  only  address  the  CPW,  and  not  the  grounded  CPW,  the  sap¬ 
phire  layer  in  our  particular  device  is  far  too  thick  to  contribute  significantly  to  the  geometric 
capacitance.  Additionally,  there  is  no  geometric  contribution  to  the  resistivity  of  the  trans¬ 
mission  line;  the  loss  induced  in  the  superconductor  is  purely  due  to  the  two-fluid  model. 
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and  will  be  taken  into  account  via  the  surface  impedance  in  a  later  chapter,  as  will  the 
kinetic  inductance. 


3.2  Lumped-Element  Model  of  a  Resonator 

The  simplest  model  of  a  resonator  comes  from  its  lumped-element  representation,  which  we 
treat  in  this  section.  It  provides  insight  to  how  the  resonator  couples  with  the  outside  world, 
in  addition  to  giving  basic  formulations  for  Q  and  oj.  In  this  section,  we  also  introduce  the 
idea  of  the  loaded  quality  factor  Ql,  which  is  the  measurable  figure  of  merit  in  this  thesis, 
as  it  includes  both  the  losses  intrinsic  to  the  resonator  and  the  losses  it  experiences  through 
coupling  with  the  outside  world. 

3.2.1  Quality  Factor  of  a  Lumped  Element  Resonator 


I 


Z 

Figure  3-5:  Parallel  RLC  circuit  with  sunusoidal  voltage  source. 


In  order  to  understand  fully  the  dynamics  of  the  superconducting  cavity  described  in  the 
introduction,  first  we  must  address  the  properties  of  a  basic  resonator.  In  its  ideal  lumped 
element  form,  a  resonator  can  be  constructed  out  of  only  two  elements:  an  inductor  and  a 
capacitor.  Energy  is  traded  back  and  forth  between  the  two  components  of  the  system,  and 
a  probe  at  either  of  the  two  nodes  of  the  circuit  would  show  that  current  and  voltage  vary 
sinusoidally  with  time.  Such  a  resonator  is  lossless,  since  its  ideal  components  have  zero 
dissipation,  but  as  we  will  be  dealing  with  a  specific  device  we  will  skip  this  treatment  and 
go  directly  to  the  relevant  circuit:  the  parallel  RLC  resonator,  as  shown  in  Figure  3-5.  The 
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reason  the  parallel  RLC  resonator  is  appropriate  is  because  the  modes  supported  (nA/2) 
match  those  of  our  device  input. 

By  analyzing  this  form  of  resonator,  we  can  develop  an  intuition  for  how  the  TL  resonator 
will  behave.  Nominally,  we  would  like  to  store  as  much  energy  in  the  cavity  as  possible  to 
generate  a  high  electric  field.  The  voltage  source  will  apply  power  to  the  circuit,  which  will 
increase  until  it  reaches  an  equilibrium  wherein  the  energy  dissipated  each  second  is  equal 
to  that  of  the  power  delivered.  To  quantify  this,  it  is  helpful  to  define  Q,  the  quality  factor 
as 


(average  energy  stored) 
(energy  loss/second) 


W„,  +  TTe 

CO — - - 

-'loss 


(3.8) 


Where  Wm  and  We  are  the  stored  magnetic  energy  (in  the  inductor)  and  stored  electrical 
energy  (in  the  capacitor),  respectively.  Since  we  wish  to  have  high  amounts  of  energy 
stored  in  the  resonator  with  low  loss,  we  desire  a  high  quality  factor.  In  order  to  examine 
the  components  of  Q  more  closely,  we  can  look  at  Z\n,  the  input  impedance  of  the  lumped- 
element  resonator,  where 


and  for  a  resonator  of  length  I,  in  terms  of  the  transmission  line  parameters 

r  _  2^1^ 


(3.9) 


(3.10) 


Now  if  we  apply  a  sinusoidal  voltage  source  across  the  two  nodes  of  the  resonator,  the 
complex  power  of  the  circuit  becomes 


Pin- 


1  1 
r^JJl 


(3.11) 
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and  from  there  it  is  straightforward  to  derive  the  components  of  Q[21].  In  particular, 


loss  2  R 


(3.12) 


(3.13) 


and 


iill! 

4a;2L’ 


finally  yielding 


Q 


2W^ 


P\c 


R 

uqL 


—  ujqRCj 


(3.14) 


(3.15) 


where  uq  —  Ify/LC  is  the  resonant  frequency  of  the  system.  This  formulation  of  Q  shows 
that  for  a  parallel  RLC  circuit,  the  quality  factor  is  directly  proportional  to  R,  which  is  the 
expected  result  since  for  P  — )•  oo  the  circuit  becomes  an  ideal  LC  circuit  with  no  dissipative 
element  and  Q  ->  oo.  It  is  also  important  to  note  that  this  quality  factor  only  describes  a 
resonant  circuit  completely  disconnected  from  the  outside  world-that  is,  when  there  is  no 
load  on  it.  When  the  circuit  is  “loaded”  with  non-ideal  external  circuitry,  invariably  the 
overall  quality  factor  will  lower  due  to  extra  loss  in  the  loading  components.  This  can  be 
seen  in  the  relation  for  the  loaded  Q,  expressed  as 


1  _  _J_  1 

Ql  Qext  Qint 

3.2.2  Coupling  to  the  Resonator 


(3.16) 


11 - ^ - II 


Figure  3-6:  Parallel  RLC  circuit  with  input  coupling  circuit  (C^,  Pl) 
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Figure  3-7:  The  same  parallel  RLC  circuit  with  norton  equivalent  input  coupling  (C*,  R*) 

In  order  to  perform  measurements  on  a  resonator  circuit,  we  must  somehow  couple  our 
measurement  input  and  output  to  the  device.  In  the  case  of  the  device  we  are  examining,  this 
coupling  can  be  described  by  a  pair  of  input  capacitors  attached  to  ground  by  load  resistors 
as  shown  in  Figure  3-6.  Though  only  a  single  resistor  and  capacitor  have  been  added  to 
each  side,  there  is  no  longer  a  trivial  way  to  extract  the  quality  factor  of  the  loaded  system. 
Fortunately,  Goppl  et  al.  have  broken  the  problem  down  into  a  much  simpler  system  that 
has  shown  good  agreement  with  experimental  results  for  device  geometries  similar  to  the 
one  described  in  this  report  [9].  By  performing  a  Norton  equivalency  of  the  coupling,  the 
series  resistor  Rl  and  capacitor  can  be  equivalently  viewed  as  a  parallel  capacitor  C* 
and  R*  as  seen  in  Figure  3-7,  specifically 


R* 


l+u;lC^,Rl 


(3.17) 


C* 


C. 


(3.18) 


where  Un  —  nujo  —  ^IVLnC  is  the  angular  frequency  of  the  nth  mode  of  resonance.  Now, 
by  summing  the  parallel  capacitors  and  parallel  resistors,  the  circuit  can  be  reformulated  as 
a  single  RLC  parallel  resonator  whose  quality  factor  represents  that  of  the  loaded  system. 


Ql. 


Ql 


,  C  +  2C* 
‘^^l/R  +  2/R* 


~  (jJ 
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c 

1/R  +  2/R* 


(3.19) 


Note  that  C*  can  be  neglected  in  the  numerator  only  because  of  the  particular  coupling 
setup  used  in  this  device  (i.e.  for  our  device  C*  C).  In  addition  to  approximating  Qlj  this 
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model  shows  that  increasing  the  coupling  capacitance  decreases  the  resonance  frequency, 
because 


^Ln{C  +  2C*y 


(3.20) 


Now  it  is  possible  to  derive  the  relationship  between  the  total  quality  factor  Qi,,  the  intrinsic 
quality  factor  Qi^t,  and  the  external  quality  factor  Qext-  In  the  literature  [21],  we  see  that  for 
a  parallel  resonant  circuit  attached  to  a  single  external  load  resistor  Ri^,  the  external  quality 
factor  can  be  shown  to  be  Qext  =  R\./{^qL)-  By  interpreting  our  external  coupling  circuitry 
as  a  single  resistor  (remember  that  C  +  2C*  «  C  so  we  can  neglect  the  external  capacitors) 
the  corresponding  load  resistor  becomes  jRl  =  R*  l‘^i  nnd  we  come  to  the  conclusion  that 


(^int  —  ^nR^  1  and 


(3.21) 


Qext  =  (3.22) 

Thus  there  are  two  regimes  of  Ql  in  which  the  coupled  oscillator  operates,  the  over¬ 
coupled  and  undercoupled.  In  the  former,  Qext  <  Qint  and  the  loaded  quality  factor  is 
governed  by  Qext  and  so  can  be  well-approximated  by  Ql  ~  C'/2a;„i?LC'^,  or  more  simply 
the  loaded  quality  factor  is  dominated  by  the  coupling  capacitance  Ql  oc  However,  at 
very  small  coupling  capacitances  Qint  <  Qext,  the  undercoupled  regime  dominates,  and  Ql 
saturates  as  can  be  seen  in  Figure  3-8.  One  interpretation  of  these  two  regimes  is  that  in  the 
overcoupled  scenario,  most  of  the  energy  the  cavity  loses  is  being  sent  out  to  the  apparatus 
that  connects  to  it.  In  the  undercoupled  scenario,  the  energy  loss  is  dominated  by  dissipa¬ 
tion  intrinsic  to  the  resonator,  such  as  radiation  or  dielectric  losses.  Experimentally,  this 
means  that  by  continually  lowering  the  coupling  capacitance,  once  we  observe  the  loaded 
quality  factor  saturating,  the  Qext  term  has  becomes  negligible  and  Qint~a  parameter  which 
is  intrinsic  to  the  resonator  and  which  need  only  be  found  once-is  determined.  Additionally, 
a  direct  relationship  is  formed  between  the  Ql  and  our  original  coupling  capacitance, 
by  Eq.  (3.16),  (3.17),  and  (3.22)  and  so  by  either  simulating  or  measuring  Ql,  the  other 
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parameter  may  be  determined. 


Figure  3-8:  Input  coupling  versus  quality  factor  in  A/2  resonator.  Blue  line  shows  the 
theoretical  relationship  between  coupling  capacitance  and  quality  factor.  Green  dashed 
lines  are  over  (sloped  line)  and  undercoupled  (horizontal  line)  regimes 


3.3  The  Superconducting  Resonator 

Up  until  this  point  we  have  treated  the  resonator  conductors  as  though  they  were  purely 
classical,  ignoring  the  effects  of  superconductivity.  In  this  section  we  explore  the  effects 
that  superconductivity  adds  to  our  model-namely  kinetic  inductance.  We  also  introduce 
the  concept  of  surface  impedance,  and  examine  the  temperature  dependence  of  Q  in  the 
superconducting  regime. 
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3.3.1  Taking  Superconductivity  into  Account 

Having  accounted  for  factors  that  affect  the  resonator  when  viewed  as  normal  metal,  we  now 
examine  the  influences  on  the  system  that  superconductivity  will  have  on  the  circuit  param¬ 
eters.  (For  simplicity,  we  assume  that  there  is  zero  magnetic  field  penetrating  the  device 
when  it  is  brought  below  the  superconducting  transition  temperature  Tc).  We  can  begin  by 
modeling  the  superconductor  as  a  lossless  conductor  and  drop  the  Ri  term  from  the  trans¬ 
mission  line  model.  However,  that  would  neglect  the  implications  of  the  penetration  depth 
A  of  the  superconductor  and  its  interaction  with  the  two-fluid  model  of  superconductivity. 
The  magnetic  field  penetration  into  the  surface  of  the  niobium  excites  both  the  Cooper  pairs 
in  the  superconducting  channel  and  the  unpaired  electrons  in  the  normal  channel.  There 
are  two  consequences  of  the  two-fluid  model  that  affect  the  superconducting  transmission 
line  model.  First,  a  kinetic  inductance  term  L^i  must  be  added  to  the  existing  geometric 
inductance  to  account  for  the  ‘inertial’  movement  of  the  oscillating  Cooper  pairs  in  the 
superconducting  channel  so  that 


Li  =  L^  +  (3.23) 

Secondly,  the  resistive  term  R\  must  be  reinstated  to  account  for  losses  incurred  by 
scattering  in  the  normal  channel.  Since  both  terms  rely  on  the  two-fluid  model,  they  will 
depend  fundamentally  on  the  ratio  of  Cooper  pairs  to  unpaired  electrons,  which  in  turn 
is  determined  by  the  ratio  T/Tc-  While  it  may  seem  that  these  terms  should  be  found 
empirically  by  varying  the  temperature  (they  can  be),  we  can  instead  relate  them  to  the 
surface  impedance  Zs  of  the  superconductor. 

3.3.2  The  surface  impedance  of  a  superconductor 

Due  to  the  magnetic  penetration  into  the  surface,  a  superconductor  can  be  thought  of  as  a 
perfect  conductor  with  a  surface  impedance 

Zs  =  R^+ju;L,.  (3.24) 
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For  simple  transmission  lines  such  as  the  parallel-plate  TL,  and  for  specific  types  of  geometry 
and  material  combinations  (such  as  thick  Al[15])  complete  analytical  solutions  exist  for  Zg. 
However,  for  the  general  case  of  a  CPW  made  of  arbitrary  materials,  the  governing  equations 
are  extremely  complex  integrals  and  must  be  solved  numerically.  To  give  an  idea  of  the 
involvement  of  the  solutions  at  low  enough  temperatures  or  at  high  enough  frequencies. 
Ohm’s  law  becomes  non-local[6].  What  this  means  is  that  the  current  density  J  at  a  point 
r  is  not  only  directly  affected  by  the  electric  field  at  that  point,  but  is  also  influenced  by 
a  complex  weighted  average  of  £  in  a  volume  around  f.  The  full  electrodynamic  analysis 
is  available  in  Gao  [8],  but  we  will  presently  skip  the  (lengthy)  details  and  instead  focus  on 
the  effects. 

We  have  coded  a  numerical  solver  for  the  surface  impedance  equations  based  on  the 
Mattis-Bardeen  kernel  formulation  originating  in  Popel[20],  and,  using  only  material  pa¬ 
rameters,  are  able  to  determine  Zg  based  the  independent  variables  of  u  and  T.  The  code 
may  be  viewed  in  Appendix  B.  The  parameters  required  are  the  London  penetration  depth 
Az,o  at  zero  temperature,  Tc,  the  mean  free  path  Imfp,  and  two  of  the  following  three  pa¬ 
rameters:  the  coherence  length  ^o,  Fermi  velocity  uq,  and  zero-temperature  binding  energy 
Ao  since  they  are  directly  related  by 


These  parameters  for  Nb[20]  are  listed  in  Table  3.3.2  for  reference. 

Table  3.3.2:  Parameters  for  bulk  Nb 


Tc  [K] 

9.2 

\lo  [nm] 

33.3 

vq  [10®  m/s] 

0.28 

Co  [nm] 

39 

Imfp  [nm] 

20 

Ao  [meV] 

1.395 

(3.25) 


Once  we  have  the  material-dependent  surface  impedance,  it  is  possible  to  determine  the 
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contribution  of  to  the  total  resonator  inductance  L  because  they  are  linearly  related  by 
a  geometric  factor  ^[Ib].  Specifically, 


Lki  =  gLs.  (3-26) 

While  normally  requiring  a  numerical  solution  of  a  contour  integral,  in  the  case  of  t  <  u) 
(such  as  our  CPW  is)  the  calculation  for  g  has  been  derived  by  Collin  [7]  and  is  estimated 
to  be  accurate  to  within  10%  for  t  <  0.025w  and  k  <  0.8.  This  approximation  for  g  is 
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(3.27) 


where  k  =  a/b. 

Now  that  we  can  calculate  p,  we  can  take  the  correct  proportion  of  Lg  and  finally  arrive 
at  Lki-  In  addition,  we  can  now  calculate  another  useful  parameter,  the  ratio  of  kinetic 
inductance  to  total  inductance 


a  = 


gLs 

Lm  +  gLs 


(3.28) 


This  parameter  describes  how  sensitive  factors  that  rely  on  the  circuit’s  inductance-particularly 
the  quality  factor  and  the  resonance  frequency-will  be  to  variables  that  affect  the  kinetic 
inductance,  most  notably  the  temperature. 

The  calculation  of  a  can  be  verified  by  comparing  the  experimental  data  of  the  under¬ 
coupled  Ql  (effectively  Qint)  with  the  quality  factor  computed  from  the  numerical  surface 
impedance  for  a  range  of  temperatures  below  Tc-  Since  the  Q  and  w  of  a  parallel  RLC 
circuit  can  be  written  as  Q  =  R/loL,  we  can  write  an  equation  for  Q{T)  and  then  fit  it 
to  the  experimental  data  at  different  temperatures  with  two  parameters;  a,  and  a  limit¬ 
ing  quality  factor  of  Q(0)  that  represents  the  baseline  power  leakage  of  the  device  at  zero 
temperature[15].  It  takes  the  form 
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Figure  3-9:  Change  in  the  inverse  quality  factor  versus  temperature 


1  1  ^oLsjT)  /n  nq'v 

^  -  Q(0)  +  aUT)  ■ 

.  Similarly  for  the  shift  in  resonant  frequency, 

Suo  _  u;o(T)  —  c<;o(0)  _  a  Ls{T)  —  Ls{0)  •jq'j 

wq  ^o(O)  2  ^s(O) 

One  such  plot  derived  from  our  surface  impedance  numerical  solver  is  shown  in  Figure  3-9. 
This  figure  shows  how  the  dependence  of  the  quality  factor  on  temperature  decreases  as  the 
temperature  approaches  zero.  Additionally,  the  higher  the  kinetic  inductance  fraction  the 
more  affected  the  quality  factor  is  by  temperature. 
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3.4  Simulation  Based  on 


With  access  to  a  numerical  solution  of  the  surface  impedance,  a  straightforward  next  step 
that  does  not  require  the  use  of  simulation  software  is  to  use  the  ABCD  method  to  analyze 
the  resonator.  Using  this  formulation,  it  is  possible  to  gain  an  understanding  of  the  current 
and  voltage  fields  in  any  two-port  transmission  line. 


3.4.1  The  ABCD  Matrix 

The  transmission  (ABCD)  matrix  is  a  simple  method  by  which  to  calculate  the  overall 
impedance  properties  of  circuit  networks,  or  in  our  case,  non-uniform  transmission  lines. 
Taking  a  section  of  a  transmission  line  Az  long,  consider  the  set  of  conductors  on  either 
edge  a  “port”  It  is  then  clear  from  Figure  3-1  that  our  section  of  transmission  line  can  be 
thought  of  as  a  two-port  network.  Using  the  ABCD  method,  the  voltage  and  current  of  one 
end  can  then  be  related  to  the  those  of  the  other  end  by  a  matrix  of  four  elements  in  the 
form 


Ui  A  B  Vi 

h  C  D  h 

where  the  elements  A,  B,  C,  and  D  can  be  calculated  from  the  circuit  model  Figure  3-1. 
The  result  is  a  straightforward  set  of  equations  that  concretely  relate  the  two  sides  of  a 
section  of  transmission  line,  and  can  be  cascaded  through  multiple  sections.  Assuming  a 
uniform  transmission  line,  in  the  limit  of  Az  0  the  original  impedance  of  the  line  would 
be  retrieved  from  the  cascaded  matrix,  along  with  the  voltage  and  current  characteristics 
at  any  point.  There  is  a  limiting  factor  that  the  transmission  line  cannot  branch  in  this 
simple  two-port  formulation,  but  the  advantage  lies  in  the  simplicity  of  the  formulation. 


(3.31) 


3.4.2  Sample  Application  for  ABCD  Superconductor  Simulation 

One  application  of  the  ABCD  method  is  to  use  a  superconducting  resonator  as  a  location- 
dependent  temperature  sensor  in  a  cryogenic  environment.  As  can  be  seen  in  Figure  3-10, 
in  a  A/2  capacitatively-coupled  resonator  such  as  ours,  in  its  first  resonant  mode  the  current 
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Figure  3-10;  Simulation  of  current  density  at  the  first  resonant  frequency  of  a  capacitatively- 
coupled  A/2  resonator 


is  most  excited  in  the  TL’s  center  (the  node  of  V^— and  therefore  antinode  of  I  in  Figure  3 
2).  At  this  point,  the  inductance  and  resistance  is  maximally  related  to  the  temperature-a 
warmer  temperature  with  a  ratio  of  T/Tc  closer  to  unity  will  have  fewer  Cooper  pairs 
providing  kinetic  inductance,  and  more  unpaired  electrons  producing  scattering  and  loss. 
If  the  resonator  were  bent  into  a  U  shape,  the  center  portion  at  the  bottom  of  the  U 
could  be  used  as  a  geometrically  static  temperature  probe.  Though  it  is  intuitive  that  the 
resonator’s  impedance— and  therefore  resonant  frequency  and  quality  factor— would  change, 
it  is  not  obvious  how  to  quantify  that  change.  It  could,  however,  be  simulated  using  the 
ABCD  method.  To  do  so,  one  could  model  a  smooth  temperature  change  over  the  length 
of  the  resonator  and  include  those  variable  values  in  the  surface  impedance  calculation  for 
each  subsection  of  the  transmission  line. 


36 


Chapter  4 


Early  Microwave  Tests  and 
Simulation 


Up  until  this  point  we  have  explored  the  physical  processes  by  which  we  fabricated  the 
resonators,  and  gained  a  theoretical  understanding  of  their  behavior.  In  this  chapter,  we 
examine  some  of  the  early  steps  undertaken  at  the  beginning  of  the  project.  Before  building 
microresonators,  we  first  fabricated  and  characterized  PCB  resonators  for  insight.  We  also 
performed  simulation  on  these  resonators  which  gave  us  a  good  understanding  of  how  the 
fields  of  the  resonators  worked. 

4.1  PCB  Resonators 

In  order  to  confirm  our  understanding  of  how  the  transmission-line  resonators  functioned, 
we  designed  a  set  of  5x10  cm  resonators  and  had  them  patterned  in  copper  on  PCB  for 
testing  at  room  temperature.  We  designed  nine  resonators  in  total,  with  varying  coupling 
capacitances,  meanders,  and  overall  lengths.  In  the  end,  one  of  the  designs  was  of  particular 
interest  to  us:  a  short  straightline  (no  meander)  resonator  with  plated  vias,  designed  to  have 
its  first  resonance  at  6  GHz.  It  is  shown  in  Figure  4-1 

We  had  the  resonators  fabricated  on  0.76  mm  Rogers  4350B^,  a  material  which  is  de¬ 
signed  for  microwave  applications  and  has  a  dielectric  constant  of  3.48.  Both  the  patterned 
^Prom  Hughes  Circuits 
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Figure  4-1;  PCB  resonator  schematic  drawing 

layer  on  top  and  the  bottom  ground  layer  were  “2  oz”  copper,  or  71  pm  thick.  The  res¬ 
onators  were  characterized  by  the  HP  8722C  network  analyzer  (detailed  in  chapter  5),  which 
connected  to  the  resonators  by  means  of  DC-50  GHz  end-launch  connectors.  The  vias  were 
drilled  into  the  resonator  board  and  plated  with  copper  on  the  recommendation  of  the 
end-launch  manufacturer,  Southwest  Microwave.  They  stated  that  the  lateral  walls  formed 
by  the  vias  on  either  side  of  the  strip  transmission  line  acted  to  prevent  additional  waveg¬ 
uide  modes  from  being  excited;  the  vias  created  what  is  called  a  “channelized  coplanar 
waveguide”  structure. 

On  testing,  we  saw  a  first  resonance  at  5.39  GHz  and  a  second  resonance  at  10.75  GHz, 
well  off  from  the  expected  6  and  12  GHz,  but  reasonable  if  the  dielectric  constants  used 
were  inaccurate.  The  quality  factor  values  were  98.5  and  102.1  for  the  first  and  second 
resonances  respectively.  At  this  point  we  needed  a  method  to  determine  if  our  results  were 
accurate.  We  decided  to  use  microwave  simulation  software  HFSS^  and  Sonnet  to  model 
the  PCB  resonator  and  see  if  matched  our  experimental  values. 

^By  Ansoft  corporation 
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4.2  Simulating  the  PCB  Resonators 

Once  we  had  the  fundamentals  of  the  general  resonator  geometry  and  testing  done,  we 
began  the  task  of  trying  to  reproduce  the  acquired  data  in  simulation  software.  We  used 
two  commercial  microwave  simulation  tools,  HFSS  and  Sonnet.  Both  operate  using  the 
finite-element  method  and  so  are  incapable  of  reproducing  transient  phenomenon,  but  that 
functionality  was  not  important  for  our  applications  because  we  were  focused  on  the  steady- 
state  phenomena  of  the  resonators. 

Though  the  two  programs  use  the  same  fundamental  technique  for  solving  fields,  they 
differ  greatly  in  that  HFSS  is  by  default  a  3D  structural  simulator,  and  Sonnet  works  by 
“stacking”  2D  layers  upon  one  another.  In  the  end,  we  chose  to  use  primarily  HFSS  because 
there  was  no  good  analog  for  building  plated  vias  between  layers  in  Sonnet  to  create  the 
channelized  CPW.  In  HFSS,  however,  we  were  able  to  define  an  artificial  perfect  electric 
conductor  boundary  that  would  function  much  the  same  way  as  the  the  plated  walls  of  the 
via  would. 

4.2.1  Building  the  resonator  model 

The  model  for  the  resonator  was  straightforward  to  build  in  HFSS,  since  the  software  comes 
with  a  set  of  predefined  dielectrics  and  conductors.  We  exported  the  layers  from  the  original 
design  then  extruded  them  to  the  appropriate  thicknesses  to  recreate  the  device  in  3D.  A 
zoomed-out  view  of  the  model  is  shown  in  Figure  4-2.  However,  defining  the  excitation  for 
the  resonator  was  not  so  simple.  Initially  we  excited  the  device  using  a  wave  port,  a  2D 
rectangle  that  imitates  the  end  of  an  infinitely  long  rectangular  transmission  line.  Later 
tests  revealed,  however,  that  changing  the  size  and  aspect  ratio  of  this  wave  port  altered 
the  resulting  scattering  parameter  profile  significantly,  as  can  be  seen  in  Figure  4-3.  The 
overarching  result  was  that  any  profile  could  have  been  fit  by  altering  the  shape  of  the 
rectangular  wave  port;  there  was  no  sensible  pattern  of  convergence. 

To  address  the  variability,  we  found  that  using  a  “lumped  port”  excitation  gave  a  much 
more  consistent  frequency  profile.  The  HFSS  documentation  also  mentions  that  the  lumped 
port  is  more  accurate  when  either  trying  to  create  a  local  excitation  on  the  PCB,  or  making 
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Figure  4-2:  HFSS  model  of  PCB  resonator 
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Figure  4-3:  Transmission  spectra  for  a  resonator  with  varying  waveport  sizes 
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electrical  contact  to  the  device  with  a  ground-signal-ground  probe,  which  is  how  the  end- 
launch  connectors  functioned.  Had  we  instead  modeled  the  complete  system,  including  the 
end-launch  connectors,  the  wave  port  would  have  been  more  appropriate.  The  edges  of  the 
2D  wave  port  sheet  function  as  a  reference  ground,  which  can  be  aligned  with  the  coax 
shield,  and  the  inner  area  would  excite  the  center  conductor. 

While  the  wave  port  excites  the  line  as  a  waveguide,  the  lumped  port  acts  as  an  excitation 
bridge  between  the  ground  plane  of  the  resonator  and  the  center  electrode.  Additionally, 
instead  of  defining  a  shape  and  allowing  the  program  to  calculate  the  entrance  impedance- 
such  as  is  the  case  with  a  wave  port-the  impedance  must  be  predetermined.  We  set  this 
impedance  value  to  50  D  and  placed  the  lumped  port  between  the  ground  plane  and  the 
input  and  output  transmission  lines,  as  shown  in  Figure  4-4. 


Figure  4-4:  Closeup  of  the  the  lumped  port  excitation  sheet  within  HFSS. 

Due  to  the  nature  of  the  meshing  in  HFSS,  we  also  defined  an  “airbox”  around  the  model 
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so  that  it  would  solve  for  the  field  above  the  resonator.  The  airbox  extended  at  least  A/4  in 
every  direction  from  the  center  electrode,  and  all  of  its  faces  (save  the  one  coinciding  with 
the  dielectric)  were  made  into  radiation  boundaries,  in  order  to  minimize  any  reflections 
stemming  from  the  mesh  ending  abruptly. 

4.3  Comparison  of  Simulation  to  Experimental  Results 

Once  we  implemented  the  lumped  port  in  our  HFSS  model,  the  results  of  the  simulation 
matched  the  PCB  resonator  extremely  well.  A  comparison  between  the  experimental  and 
simulated  data  is  shown  in  Figure  4-5  .  The  resonances  are  clearly  visible  in  the  simulation, 
their  frequencies  matching  the  experimental  data  to  within  10%  and  their  quality  factors 
to  within  5%.  The  remainder  of  the  spectrum  away  from  resonance  in  both  cases  is  -20  dB 
lower  than  the  resonant  peaks,  or  less  than  1%  of  the  maximum  output  power. 


Figure  4-5;  Comparison  of  measurement  (blue)  and  simulation  (red)  of  the  PCB  resonator. 
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4.3.1  Simulation  field  distribution 


With  the  simulation  confirming  our  experimental  results,  the  next  step  was  to  examine  the 
microwave  field  distribution  at  resonance.  We  looked  at  two  slices  through  intersecting  the 
symmetrical  origin  of  the  resonator:  one  perpendicular,  and  one  parallel  with  the  center  of 
the  resonator  electrode.  As  can  be  seen  in  Figure  4-6,  when  exciting  the  resonator  into  its 
second  resonance,  the  parallel  slice  shows  the  sinusoidal  variation  in  magnitude  one  would 
expect  in  the  electric  field.  The  anti-nodes  are  clearly  visible  in  the  center  and  at  either 
coupling  capacitor,  forming  one  period  of  a  cosine,  just  as  transmission  line  theory  would 
suggest. 


A 


Figure  4-6:  PCB  electric  field  plotted  in  HFSS  at  second  resonance 


In  the  perpendicular  slice,  the  symmetry  of  the  system  produces  a  weak  z-directed  field 
in  the  center,  and  higher  fields  along  the  shorter  distance  to  ground:  between  the  electrodes 
of  the  top  layer  of  the  PCB.  Since  the  Rogers  dielectric  below  the  CPW  is  0.76  mm  thick. 
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the  majority  of  the  field-and  hence  the  energy  in  the  mode-is  carried  between  the  coplanar 
ground  and  electrode,  which  are  only  0.28  mm  apart. 


4.4  Simulating  the  Niobium  Microresonators 

Once  the  PCB  resonators  had  been  simulated,  the  next  step  was  to  adapt  the  geometry 
of  the  model  to  fit  that  of  the  actual  superconducting  resonators  we  fabricated.  Although 
at  the  time  we  were  fabricating  meandered  resonators,  the  curves  of  the  electrodes  greatly 
increased  the  amount  of  mesh  points  in  the  model.  Instead,  we  opted  to  simulate  the 
straight-line  resonator  to  reduce  the  amount  of  time  and  memory  necessary  to  run  the 
simulation.  This  also  simplified  the  visual  inspection  of  the  electric  field  once  the  simulation 
was  complete. 

4.4.1  Adapting  the  PCB  model 

We  created  the  microresonator  model  in  the  same  way  as  that  of  the  PCB  resonator,  only 
changing  geometric  values,  materials,  and  removing  the  perfect  electric  boundaries  since 
there  was  no  equivalent  to  the  PCB  vias  in  the  sapphire-niobium  fabrication.  Unlike  the 
PCB  resonator,  the  channelization  of  the  niobium  microresonator  CPW  was  unecessary 
because  the  sapphire  dielectric  was  so  thick  compared  to  the  resonator  gap  spacing  that 
no  unwanted  waveguide  mode  existed  in  the  frequency  range  of  interest.  The  dielectric 
was  changed  from  PCB  material  to  sapphire,  with  a  dielectric  constant  of  11.5  correspond¬ 
ing  to  the  C  axis  of  the  sapphire.  Additionally,  the  conductor  was  changed  to  a  lossless 
perfect  conductor.  By  just  using  the  perfect  conductor,  this  model  neglected  to  include 
kinetic  inductance  of  the  surface  impedance  for  the  sake  of  simplification.  However,  be¬ 
cause  the  320  nm  niobium  layer  is  many  times  thicker  than  the  London  penetration  depth 
of  47±5  nm[14],  the  kinetic  inductance  contribution  to  the  overall  inductance  is  small. 

The  primary  challenge  with  simulating  the  microresonators  was  refining  the  mesh  while 
remaining  under  memory  limits  of  the  modeling  computer.  It  was  particularly  difficult 
because  of  the  large  difference  in  length  scales  between  features:  the  length  of  the  entire 
device  is  over  20  mm,  but  the  smallest  feature  size  (the  gap  capacitors)  was  5  pm.  The 
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HFSS  meshing  process  puts  more  mesh  points  around  areas  of  fine  detail,  but  the  four  orders 
of  magnitude  difference  between  our  length  scales  caused  it  problems.  We  had  to  define  a 
mesh  refinement  box  localized  above  and  around  the  resonating  electrode  that  forced  fine 
meshing,  in  order  to  get  the  scattering  parameters  to  converge.  The  results  for  the  field 
just  above  the  resonator  electrode  and  gaps  at  second  resonance  can  be  seen  in  Figure  4-7. 


Figure  4-7;  Magnitude  of  the  electric  field  versus  XY  position  at  second  resonance. 

Additionally,  we  plotted  the  field  dependence  along  the  Z  direction.  As  seen  in  Figure  4- 
8  the  field  falls  off  exponentially  with  distance  from  the  resonating  electrode.  This  matched 
our  expectations,  and  served  as  a  sanity  check  for  the  model. 

4.4.2  Simulating  the  coupling  capacitances 

In  order  to  compare  how  the  physical  devices  acted  with  the  lumped-element  model,  we  first 
needed  to  simulate  the  coupling  capacitances  for  each  resonator.  Each  coupling  capacitor 
was  simulated  in  HFSS  by  modeling  the  geometry  as  shown  in  Figure  4-9,  from  1  to  10  GHz. 
Since  the  electrode  material  was  selected  to  be  a  perfect  electric  conductor,  the  admittance 
was  completely  imaginary  and  capacitative;  Y21  —  The  resulting  admittance  values 

Y21  were  then  used  to  find  the  capacitance  at  any  frequency. 


46 


b _ I _ I _ I - \ - 1 - 1 - 1 — 

0  0.01  0.02  0.03  0.04  0.05  0.06  0.07 


Distance  above  resonator  surface  (mm) 

Figure  4-8:  Log  magnitude  of  the  electric  field  versus  distance  from  middle  of  the  resonator 
electrode  at  its  second  resonance.  Note  that  as  the  field  gets  further  from  the  near-field 
effects  of  the  electrodes  it  converges  to  an  exponential  fallolT  as  expected 


Figure  4-9:  Coupling  capacitor  modeled  in  HFSS 
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4.4.3  Including  kinetic  inductance  in  the  model 

The  next  step  for  the  simulation  model  is  including  the  kinetic  inductance  portion  to  produce 
more  accurate  results.  This  can  be  accomplished  through  the  judicious  use  of  surface 
impedance  along  the  faces  of  the  box  that  define  the  resonator  center  electrode.  The  surface 
impedance  will  be  in  the  form  of  a  resistance  plus  a  reactance  [28]  [20],  and  can  be  calculated 
from  material  parameters  as  is  addressed  in  section  3.3.2. 


48 


Chapter  5 


Testing  and  Characterizing  the 
Superconducting  Resonators 


In  this  chapter  we  describe  the  testing  procedures  used  when  characterizing  the  super¬ 
conducting  resonators.  We  used  two  sets  of  testing  equipment:  we  found  niobium  material 
parameters  using  a  dipstick  with  a  four-point  probe,  and  the  we  characterized  the  resonators 
with  a  cryogenic  probe  station.  We  also  detail  the  microwave  and  readout  electronics  used 
in  conjunction  with  the  probe  station. 

5.1  Dipstick  Material  Measurements 

In  the  following,  we  describe  the  procedure  we  used  to  determine  Tc  for  our  niobium. 
Additionally,  we  describe  the  method  by  which  the  dipstick  testing  apparatus  can  determine 
the  residual  resistance  ratio  and  sheet  resistance  for  our  niobium  film.  Though  these  latter 
tests  were  never  executed  to  conclusion,  they  will  be  in  the  future  and  so  appear  here  for 
reference. 

5.1.1  Measurement  of  Tc 

In  order  to  characterize  the  superconducting  properties  of  the  material  stack,  we  ran  a  four- 
point  measurement  on  unpatterned  samples  as  they  cooled  to  liquid  helium  temperatures. 
The  sample  was  mounted  to  the  end  of  a  dipstick,  which  was  then  submerged  in  a  helium 
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bath  inside  of  a  dewar.  Four  evenly  spaced  gold  spring  pins,  arranged  in  a  line,  pressed 
against  the  niobium  layer,  and  current  was  passed  between  the  outer  two  pins,  as  shown 
in  Figure  5-1.  This  method  had  the  advantage  of  removing  any  input  resistances  from  the 
measurement [26],  depicted  as  resistors  in  Figure  5-1.  While  the  sample  cooled  from  297K  to 
4.2K,  a  measurement  of  the  voltage  between  the  two  inner  pins  was  recorded  continuously, 
and  the  resistivity  of  the  sample  could  be  calculated  by  factoring  in  the  relative  pin  distance. 


Figure  5-1:  Four  point  measurement  depicting  current  flow,  measurement  voltmeter,  and 
pin  alignment 

The  general  form  for  the  resistivity  derived  from  the  four-point  measurement  is 

P  =  Gj  (5.1) 

where  G  is  a  correction  factor  dependent  on  sample  shape,  dimensions,  and  the  arrange¬ 
ment  of  the  four  electrical  contacts.  With  the  dipstick,  we  were  measuring  an  extremely  thin 
conductor  with  the  four  electrical  contacts-spring  pins-in  a  line,  and  the  thickness  of  the 
niobium,  t,  is  significantly  smaller  than  the  spacing  between  spring  pins,  s.  The  resulting 


50 


G  is  given  [25]  by 


G  =  ~tT2{-)  -  4.5324i  (5.2) 

In  2  s 

since  22(|)  — )■  1  as  |  0. 

There  is  an  additional  caveat,  however.  This  correction  factor  assumes  that  the  metal 
layer  is  infinitely  long  and  infinitely  wide  in- plane.  What  this  translates  to  in  practical  terms 
is  that  the  surface  under  test  had  to  have  a  diameter  of  at  least  an  order  of  magnitude  larger 
than  s  or  else  an  additional  correction  term  C  must  be  multiplied  into  G,  defined  as 


Col- 
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(5.3) 


The  easiest  way  to  get  around  this  wasto  conduct  the  four-point  measurement  in  the  center 
of  the  full  5  cm  wafer. 


5.1.2  Tc  measurement  procedure 

To  perform  the  dipstick  measurement,  first  a  sample  no  larger  than  3  cm  square  was  fixed 
to  the  sample  holder  on  the  end  of  the  dipstick  using  only  the  pressure  applied  by  the  gold 
spring-probe  pins.  With  the  sample  holder  at  25  ,  the  four  point  resistivity  measurement 
began  by  sourcing  100  /xA  of  current  across  the  outer  pins  and  logging  the  voltage  of  the 
inner  pins.  The  sample  holder  was  then  placed  into  the  dewar’s  central  chamber,  secured, 
and  the  outer  sleeve  of  the  dewar  was  filled  with  liquid  nitrogen.  The  dewar  can  be  seen  in 
Figure  5-2. 

A  probe  on  the  sample  holder  read  the  current  temperature  of  the  dipstick,  and  during 
this  process,  the  inner-pin  voltage  was  continuously  correlated  to  the  temperature.  When 
the  liquid  nitrogen  brought  the  holder  temperature  below  115  K,  we  began  adding  liquid 
helium  to  the  dewar’s  central  chamber.  We  flowed  liquid  helium  into  the  central  chamber 
until  the  probe  read  4.2  K,  at  which  point  we  continued  flowing  for  another  90  sec  to  ensure 
an  adequate  bath.  At  this  point,  the  temperature  of  the  sample  could  be  controlled  by 
raising  and  lowering  the  dipstick,  which  moved  it  through  varying  vapor  pressures  above 
the  helium  bath.  Using  this  method,  the  temperature  was  slowly  modulated  around  4  to 
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Figure  5-2:  Dewar  and  dipstick  setup  for  cooling  down  samples 
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30  K,  to  find  the  critical  current  and  ensure  that  any  lag  between  the  temperature  and 
resistivity  readings  was  negligible  versus  the  rate  of  heating  and  cooling.  Once  data  was 
done  being  collected,  the  dipstick  was  either  removed  from  the  dewar  for  sample  retrieval, 
or  left  to  warm  slowly.  Either  way,  the  sample  was  not  used  again,  as  it  had  just  endured 
cooldown  and  the  pressure  of  gold  spring-probe  pins  andwas  no  longer  a  suitable  surface  for 
fabrication. 

5.2  Device  Characterization  with  the  Probe  Station 

The  majority  of  the  testing  done  in  this  thesis  was  performed  in  the  cryogenic  probe  station. 
Although  not  as  well  shielded  as  a  complete  immersion  system,  the  probe  station  allowed  for 
rapid  testing  of  multiple  devices  on  a  single  chip.  Using  only  a  network  analyzer,  we  were 
able  to  gather  all  the  relevant  data  necessary  to  determine  quality  factors  of  the  devices 
under  test.  The  procedure  for  cooling  the  station  down  was  very  involved  and  is  described 
in  section  5.2.3. 

5.2.1  Probe  station  apparatus 

The  probe  station  we  used  for  testing  was  a  Desert  Cryogenics  TTP4,  with  a  cold  head 
capable  of  reaching  1.4  K.  The  probes  we  used  were  Picoprobe  Model  P-14-6204-04,  shown 
in  Figure  5-4.  The  system  is  pictured  in  Figure  5-3,  and  consists  of  a  cylindrical  central 
chamber  with  the  cold  head  area  in  its  center.  The  cold  head  lies  in  the  center  of  the 
chamber,  and  as  the  coldest  point  of  the  system  it  is  where  the  sample  holder  and  sample 
is  mounted.  Since  the  probe  station  is  open-cycle,  the  cold  head  is  cooled  down  with  an 
external  dewar  of  helium.  A  roughing  pump  is  also  attached  to  the  plumbing  of  the  system 
to  lower  the  pressure  beneath  the  cold  head  and  allow  the  station  to  reach  temperatures 
below  4.2  K. 

The  chamber  wasl6  cm  in  diameter  with  a  window  in  its  cap  and  bellows  protruding 
from  its  sides.  It  is  shown  open  in  Figure  5-5.  A  camera  was  positioned  above  the  system, 
pointing  down  through  the  chamber  window.  Its  lens  system  combined  viewing  and  lighting 
into  one  optical  column,  which  allowed  us  to  illuminate  and  zoom  in  on  the  areas  we  were 
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Figure  5-4;  Schematic  of  the  probes  used  in  the  TTP4  probe  station.  The  left  shows  the 
shape  of  the  body,  the  right  shows  a  closeup  of  the  probe  tip  and  its  ground-signal-ground 
setup 

probing.  Since  the  structures  we  probed  are  on  the  scale  of  microns,  the  camera  zoom  was 
necessary  to  navigate  the  probes  around  the  sample. 

Each  of  the  bellows  protruding  from  the  chamber  contained  an  arm,  with  the  probe  at 
its  end.  The  probe  was  rigidly  fixed  to  the  end  of  the  arm,  but  the  arm  was  capable  of 
movement  in  all  three  cartesian  axes,  allowing  the  probe  tip  to  contact  different  areas  of  a 
sample  fixed  to  the  cold  head.  A  stainless  steel  microwave  feedthrough  that  ran  the  length 
of  the  arm  delivers  microwave  power  from  the  outside  world  to  the  probe.  To  sink  the  heat 
carried  by  this  connection,  the  probe  body  was  connected  to  the  cold  head  by  a  thick  copper 
braid. 

Surrounding  the  cold  head  inside  the  chamber  was  the  radiation  shield,  a  stainless  steel 
radiation  barrier  with  an  infrared-absorbing  window  aligned  with  the  chamber  window  for 
visibility.  It  acted  as  a  thermal  boundary  between  the  300  K  chamber  walls  and  the  sample. 
This  radiation  shield  had  ports  in  line  with  each  bellows,  through  which  the  probes  could 
reach  the  cold  head. 


55 


^1™ 

P’:' ■ 

5.2.2  Measurement  electronics 


Through  a  GPIB  interface,  we  collected  data  about  the  system  through  three  devices.  A 
four-point  temperature  sensor  mounted  on  the  cold  head  was  connected  via  four-wire  cryo¬ 
genic  cable  ^  to  a  Lakeshore  321  temperature  controller.  An  additional  sensor,  monitoring 
the  temperature  of  the  radiation  shield  and  sample  holder,  was  connected  to  a  Lakeshore 
331  temperature  controller.  Lastly,  to  determine  the  microwave  characteristics  of  the  res¬ 
onators,  we  connected  an  HP  8722C  Network  Analyzer  to  the  SMA  connectors  on  the 
station’s  feedthrough. 

The  network  analyzer  had  two  ports,  arid  was  able  to  capture  reflection  and  transmission 
scattering  parameters  over  a  range  of  50  MHz  to  40  GHz.  Since  our  first  resonance  occurred 
around  3  GHz,  we  typically  looked  at  a  range  inclusive  of  the  first  three  resonances,  1  to 
10  GHz.  To  determine  the  quality  factor  of  high-Q  devices,  we  sampled  data  around  reso¬ 
nance  at  the  network  analyzer’s  maximum  resolution,  100  KHz.  We  collected  the  real  and 
imaginary  portions  of  the  scattering  parameters  versus  frequency.  The  data  was  collected 
using  a  Prologix  serial  GPIB-USB  Controller  connected  to  a  laptop  computer.  We  wrote 
two  layers  of  MATLAB  code,  the  first  to  interface  with  the  Prologix  interface,  and  on  top 
of  that  a  second  layer  specific  to  the  8722C,  to  facilitate  the  collection  of  the  GPIB  data. 
This  code  is  available  in  Appendix  A. 

Since  we  fabricated  A/2  resonators,  we  were  most  interested  in  collecting  the  transmis¬ 
sion  scattering  parameter  S21  with  the  network  analyzer.  Tests  showed  that  the  devices 
reflected  a  large  majority  of  the  power  incident  upon  them,  even  at  resonance,  and  so  the 
dips  that  should  be  visible  in  reflection  (Sn)  measurements  were  only  a  small  fraction  of 
the  total  power  returned,  making  them  difficult  to  analyze.  This  was  likely  due  to  mis¬ 
matched  input  impedances,  or  with  the  overall  loss  of  the  system.  Although  we  initially 
experimented  with  external  amplifiers  and  filters,  once  we  had  a  reliable  fabrication  process, 
the  transmitted  power  around  resonance  was  high  enough  to  make  external  amplification 
and  filtering  unnecessary. 

The  primary  sources  of  microwave  loss  in  the  system  were  from  the  stainless  steel 

^Prom  Lakeshore 
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feedthroughs  and  the  probes  themselves.  The  probe  insertion  loss  graph  is  shown  in  Fig¬ 
ure  5-6.  More  critically,  the  feedthroughs  have  stainless  steel  shields  and  cores  which  mini¬ 
mized  thermal  loading  of  the  probes,  but  which  are  also  highly  lossy  at  microwave  frequen¬ 
cies.  When  the  system  is  cold,  one  end  of  the  feedthrough  is  cold,  and  the  other  is  outside, 
making  it  impossible  to  isolate  the  feedthrough  for  microwave  characterization.  We  were 
able,  however,  to  place  the  probes  very  close  together  (signal  tips  within  10  microns)  on  a 
transmission  line,  essentially  shorting  probes  together.  The  result  is  shown  in  Figure  5-7, 
and  its  maximums  represent  the  total  loss  of  the  probe  station.  The  large  dips  in  the  figure 
are  due  to  the  extra  capacitance  of  the  transmission  line.  In  future  tests,  a  specific  cali¬ 
bration  structure  with  only  small  contact  pads-not  the  center  of  a  20  mm  TL-should  be 
implemented  on  the  device  surface. 


Figure  5-6:  S21  and  Sn  characteristics  of  the  probe. 


5.2.3  Testing  procedure 
Mounting  the  sample 

Finding  a  method  to  mount  the  sample  on  the  cold  head  required  a  number  of  revisions 
to  our  sample  holder.  Initially  we  began  with  a  vertically  mounted  sample,  and  fixed  a  Si 
mirror  to  the  holder  at  45°to  be  able  to  see  the  sample  with  the  camera.  This  system  was 
first  adopted  because  the  other  users  of  the  probe  station  were  using  the  vertical  mounting 
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system.  It  had  two  large  drawbacks,  however:  the  upside-down  T  shape  of  the  holder  was 
structurally  poor  at  cooling  the  sample,  and  it  was  extremely  difficult  to  align  the  Si  mirror 
with  the  camera  to  see  the  vertically-mounted  sample. 

The  first  major  revision  was  removing  the  vertical  mounting  area  and  making  the  sample 
holder  a  single  rectangular  block  of  OFHC  copper,  upon  which  the  sample  would  be  affixed 
by  conductive  silver  paint.  This  improved  the  visibility  of  the  sample  greatly  through  the 
camera,  as  we  no  longer  required  an  aligned  mirror.  Unfortunately,  for  several  tests  after 
this  change  we  struggled  with  highly  inconsistent  results.  Eventually,  we  determined  that 
the  sample  was  superconducting  only  in  certain  areas,  and  that  there  were  localized  warm 
spots  due  to  poor  and  uneven  heat  conduction  through  the  silver  paint.  We  changed  the 
type  of  silver  paint  to  SPI  Silver  Paste  Plus,  more  suitable  for  cryogenic  heat  conduction, 
and  the  testing  consistency  improved.  Eventually,  we  settled  on  creating  a  frame  that 
would  press  the  sample  down  onto  the  holder  to  be  able  to  use  indium  as  a  conduction 
layer,  since  the  silver  paste  had  a  thermal  conductivity  of  4  W/mK  versus  that  of  indium, 
81.8  W/mK.  Between  the  sample  and  sample  holder  we  put  a  layer  of  0.25  mm  indium  for 
heat  conduction.  The  holder  assembly  can  be  seen  in  Figure  5-8.  To  further  increase  the 
contact  of  the  indium,  after  the  sample  was  fixed  in  place  by  the  frame,  the  entire  holder 
setup  was  baked  at  90°  C  for  5  min. 

Cooling  down  the  probe  station 

In  order  to  cool  down  the  probe  station,  first  all  of  the  air  had  to  be  pumped  out  of  it  to 
minimize  gas  condensation.  Once  the  sample  holder  was  mounted  on  the  cold  head,  the 
turbo  was  engaged  until  the  chamber  reached  a  pressure  below  3x10“®  mbar.  Once  under 
vacuum,  we  flowed  liquid  helium  into  the  cryogen  plumbing.  Initially,  the  needle  valve 
underneath  the  station  was  closed  to  force  the  helium  to  cool  the  radiation  shield.  If  the 
cold  head  was  cooled  first,  the  remaining  gas  in  the  chamber  was  prone  to  condensing  on 
the  sample.  Once  the  radiation  shield  dropped  below  50  K,  the  needle  valve  was  opened 
slightly  and  the  helium  roughing  pump  was  turned  on  to  cool  down  the  cold  head.  The 
cooling  reached  equilibrium  when  the  cold  head  reached  1.4  K  and  the  radiation  shield 
dropped  below  30  K.  Bringing  the  probe  station  from  vented  to  cryogenic  temperatures 
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took  about  3  hours  and  10-15  L  of  liquid  helium.  The  cold  head  could  remain  stable  at 
1.4±0.1  K  at  a  cost  of  about  5  L/hour. 

Testing  the  devices 

Once  the  system  had  reached  the  desired  temperature,  the  probes  were  maneuvered  to  test 
individual  devices.  Once  aligned  with  the  device,  the  probe  signal  fingers  were  dropped  onto 
either  end  of  the  center  electrode  of  the  resonator,  and  DC  multimeter  tests  were  performed 
before  taking  any  microwave  measurements.  First,  we  measured  the  resistance  between  the 
ground  and  signal  pin  of  the  probes.  A  finite  resistance  indicated  the  presence  of  a  short. 
Next,  the  two  signal  lines  of  the  probes  were  checked  for  connectivity  to  ensure  the  center 
electrode  was  unbroken. 

Once  the  DC  tests  were  complete,  the  probes  were  moved  to  the  device  contact  pads  for 
microwave  characterization.  We  attached  the  network  analyzer  to  the  system  using  double- 
shielded  50  D  SMA  cables.  The  cable  impedances  were  calibrated  out  of  the  measurement 
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using  the  8722C  calibration  kit.  Unfortunately  due  to  the  fixed  nature  of  the  feedthroughs, 
their  loss  and  impedances  could  not  be  calibrated  against.  Typically  we  took  a  large  sweep 
with  1601  points  on  the  network  analyzer,  from  1  to  10  GHz.  At  this  resolution,  resonances 
with  a  quality  factor  of  10  to  10,000  were  plainly  visible.  The  spectrum  was  recorded  in 
real  and  imaginary  components  of  S21  for  this  range,  and  again  after  zooming  in  on  specific 
resonances.  The  results  are  discussed  in  the  following  chapter. 
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Chapter  6 


Results 


Once  we  had  fabricated  and  tested  the  devices,  we  then  needed  to  characterize  them.  This 
was  done  by  fitting  parameters  to  each  set  of  data  from  the  network  analyzer,  particularly 
Q  and  fo-  We  found  that  we  had  achieved  a  maximum  Q  of  1.2x10^  in  the  probe  station. 
Additionally,  by  simulating  the  coupling  capacitances  we  were  able  to  fit  identical  devices 
with  varying  coupling  capacitances  against  their  expected  lumped-element  behavior,  and 
found  good  agreement.  We  also  have  assembled  data  from  a  resonator  with  light  shining 
on  it,  and  examined  the  relationship  between  the  light  and  the  sample  temperature.  The 
temperature  readings  were  unfortunately  unreliable  for  much  of  the  time  during  which  this 
data  was  taken,  but  in  the  future  we  plan  to  move  to  a  cryogenic  immersion  testing  setup 
in  a  dewar,  at  which  point  we  should  have  much  more  accurate  temperature  readings  and 
control. 

6.1  Fitting  Quality  Factors 

Fitting  quality  factors  from  S21  was  performed  in  MATLAB  in  three  steps.  First,  we 
estimated  Q  using  by  setting  the  maximum  magnitude  point  as  fo,  and  then  finding  the 
halfway  point  on  either  side  of  the  peak  to  get  the  FWHM,  A/,  and  then  calculating  the 
quality  factor  using  Q  =  /o/A/.  This  estimate  of  fo  and  Q  were  then  passed  as  initial 
parameters  to  the  least-squares  nonlinear  solver  Isqnonlin,  fitting  the  equation 
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against  the  data  provided.  In  this  equation,  A,  B,  C,  fo,  and  Q  are  the  parameters  for 
the  lorentzian.  The  lorentzian  fit  given  in  the  above  equation  was  used  because  of  its  noise 
tolerance  and  excellent  results  when  compared  with  other  fitting  methods[19][4]. 

The  least-squares  solver  typically  found  a  good  fit,  but  visual  inspection  revealed  it 
often  found  a  local  minima  that  did  not  quite  align  with  the  data  points.  To  get  slightly 
better  fits,  the  resulting  parameters  from  Isqnonlin  were  then  fed  into  another  optimization 
function,  fminsearch.  The  fit  parameters  found  by  fminsearch  are  the  ones  used  in  all 
discussion  of  values  in  this  chapter. 


6.2  High-Q  Results 

Using  the  testing  setup  detailed  in  Chapter  5,  we  successfully  tested  resonators  with  a 
maximum  Q  of  1.2x10“^.  The  results  for  the  first  resonance  around  3  GHz  are  shown  in 
Figure  6-1.  Since  the  network  analyzer  had  poor  resolution  relative  to  the  bandwidth,  there 
are  few  data  points  around  the  peak,  but  the  sides  of  the  peak  were  detailed  enough  to 
generate  a  good  fit. 

The  same  resonator  is  shown  in  Figure  6-2,  demonstrating  the  circle  fit  technique  on 
the  real-imaginary  axes.  It  serves  to  highlight  how  even  a  small  number  of  points  can  be 
fitted  to  define  the  shape  of  the  entire  transmission  spectrum. 

This  resonator  had  a  center  strip  width  of  7pm  and  gaps  of  4pm,  and  the  quality  factors 
of  it  were  significantly  higher  than  those  of  the  the  resonators  with  larger  center  strip 
widths.  We  have  not  yet  gathered  data  on  center  strip  width  versus  quality  factor,  but  the 
literature[27]  suggests  that  smaller  center  strips  yield  higher  quality  factors. 


6.3  Q  vs  Ck 

By  calculating  the  distributed  element  values  as  described  in  section  3.2  and  simulating 
the  coupling  capacitances,  we  were  able  to  plot  the  expected  relationship  between  coupling 
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Figure  6-1;  First  resonance  of  the  highest-Q  resonator  we  tested.  Experimental  data  points 
are  in  blue,  the  lorentzian  fit  is  in  red  with  parameters  fo  —  3.235  GHz,  Q  =  12,630. 


capacitance  and  Q,  as  discussed  in  Chapter  3  and  compare  them  with  our  experimental 
results.  The  results  are  shown  in  Figure  6-3.  They  show  good  agreement  with  the  shape  of 
the  expected  results,  but  some  of  the  experimental  values  of  Q  were  off  by  a  factor  of  2  or 
more.  This  could  be  due  to  a  number  of  reasons,  foremost  among  which  is  the  uncertainty 
in  temperature.  At  the  temperatures  we  were  working  with,  T/Tc  0.5,  and  so  small 
variations  of  sample  temperature  translate  into  large  variations  of  Q.  Other  factors,  such  as 
altered  input  impedance  Rl  or  different  lumped-element  values  of  L  or  C  are  unlikely,  since 
they  would  not  vary  so  much  from  resonator  without  significant  fabrication  issues  which 
were  not  observed. 

The  CPW  calculations  gave  geometric  inductance  and  capacitance  per  unit  length  of 
=  8.16  nH/m  and  =  3.36  pF/m.  The  resonators  were  20.2  mm  long,  yielding  a 
total  L  of  1.64  nH  (for  the  first  resonance)  and  C  of  1.39  pF.  Examining  the  highest  values 
of  Q,  we  assume  Qint  to  be  5400,  and  from  Eq.  (3.21)  we  arrived  at  an  R  value  of  186  kfi. 

6.4  Dipstick  Cooldown  Results 

The  resistivity  measurements  during  cooldown  for  Nb  and  NbN  are  shown  in  Figure  6-4 
and  Figure  6-5  respectively.  The  graphs  clearly  delineate  the  transition  between  the  normal 
conducting  state  and  superconducting  state,  showing  us  the  critical  temperature  Tc  for 
each  material.  The  residual  resistance  ratio  was  not  found,  however.  This  is  because  an 
automatic  voltage  readout  was  unavailable  at  the  time  of  experiment  and  the  temperature 
was  not  logged  during  the  early  stages  of  cooldown.  Additionally,  the  graphs  are  in  terms  of 
the  centerpin  voltage  rather  than  the  calculated  resistivity  p;  at  the  time  of  the  experiment  it 
was  not  clear  that  the  geometry  and  dimensions  of  the  sample  needed  to  taken  into  account, 
so  arbitrary  shaped  samples  were  used  that  are  too  small  to  be  negligible  compared  to  the 
pin  spacing  s. 

6.5  Superconducting  Transmission  Line  Measurements 

As  a  basis  for  testing  in  the  probe  station,  next  to  every  set  of  resonators  we  fabricated 
we  also  put  a  transmission  line  with  the  same  geometry  (without  the  coupling  capacitors) 


67 


for  reference  purposes.  Just  like  the  resonators,  they  were  designed  with  a  characteristic 
impedance  of  Zq  =  50  fi.  The  results  are  shown  in  Figure  6-6,  and  they  indicate  the  system 
is  behaving  as  expected.  The  slope  of  the  magnitude  is  constant  with  some  fluctuation  due 
to  the  uncalibrated  impedance  mismatch  of  the  system  causing  voltage  standing  waves. 


Figure  6-6:  S21  of  a  superconducting  transmission  line.  It  has  a  length  of  20  mm.  Its 
centerstrip  has  a  gap  of  10  pm,  and  width  of  18  pm 

6.6  Results  with  Varying  Incident  Light  intensities 

To  investigate  the  effects  of  light  incident  on  the  devices,  we  took  the  variable  light  source 
from  the  camera  and  directed  it  through  the  chamber  window.  We  then  measured  the 
light  reflected  back  out  of  the  window  with  a  photodiode,  taking  that  to  be  approximately 
proportional  to  the  intensity  incident  on  the  devices.  The  results  are  shown  in  Figure  6-7. 
We  also  graphed  the  change  in  Q  and  fo  with  the  light  intensity,  shown  in  Figure  6-8  and 
Figure  6-9. 
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Figure  6-7:  Second  resonances  with  varying  intensities  of  light  incident  on  the  device 
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Figure  6-8:  Quality  factor  versus  intensity  of  light  incident  on  the  device 
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Figure  6-9:  Resonant  frequency  versus  intensity  of  light  incident  on  the  device 
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The  datapoints  in  both  graphs  have  the  same  relationship  versus  light  intensity  as  they 
do  versus  temperature,  which  led  us  to  believe  that  the  two  were  proportional.  Although  we 
did  not  have  accurate  temperature  readings  for  this  data,  we  show  this  proportionality  in 
Figure  6-10,  where  we  have  fitted  this  light-temperature  relationship  for  fo  by  hand.  Using 
computed  values  of  Lm,  Lg,  and  p,  we  determined  a  and  then  selected  the  temperature 
range  over  which  the  data  fit  Eq.  (3.30)  best,  in  this  case  3  to  5  K. 


Figure  6-10;  Change  in  resonant  frequency  /o  versus  temperature  change  caused  by  light 
incident  on  a  niobium  resonator.  In  red  is  the  calculated  A/o,  while  the  blue  circles  are 
data  points 
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Appendix  A 


MATLAB  8722C  GPIB  Interface 
Code 


%  File:  get_8722c_sweep.m 

7„  Creation  Date :  28-Sep-2009 

7,  Modification  Date:  Ol-Nov-2009 
7o  Author:  Adam  McCaughan  (amcc@mit.edu) 

7. 

7.  Grabs  whatever  data  is  on  the  screen  currently  by  switching  it  to  single 
%  mode,  collecting  the  data,  then  returning  it  to  continuous 


function  [F,M]  =  get_8722c_sweep(sl) 


%  Since  the  8722C  is  on  GPIB  address  7 
serial_write(sl,  ’++addr  7’) 


fstart  =  strread(serial_write(sl ,  ’STAR?;’)); 
fspan  =  strread(serial_write(sl,  ’SPAN?;’)); 
fstop  =  strread(serial_write(sl ,  ’STOP?;’)); 
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numpts  =  strread(serial_write(sl ,  ’P0IN?;O); 
f_inc  =  f span/ (numpts- 1) ; 

F  =  fstart+((l:numpts)-l) .*f_inc; 

F  = 

serial_write (si ,  ’ OPC? ; SING; ’ ) ; 

serial_wait(sl,300) ; 

serial_write(sl,  ’F0RM4;’); 

serial_write(sl ,  ’ OUTPFORM ; ’ ) ; 

pause  (2);  °L  pause  while  it  readies  output 

M  =  serial_readn(sl ,  numpts); 

M  =  slist2nums(M,  ’‘/»n,°/onO; 

M  =  M(: ,1); 


/O  /0/0/9/9/8  /o  /o 


0/  0/  0/  0/  0/  0/  0/  0/ 


%  File;  run_8722c_sweep .m 
%  Date:  28-Sep-2009 

y.  Author:  amcc@mit.edu 

I 

y.  Takes  in  a  frequency  start  (GHz)  and  stop,  number  of  points  to  sample 
y,  between  those  two,  and  returns  the  frequency  &  magnitude  data  from  them 
y.  Needs  si  the  serial  object 

y. 

y.  Returns  magnitude  data  in  the  format  [Magi  Mag2  .  .  .  MagN]  where  N  =  navgs 


function  [F,M]  =  run_8722c_sweep(sl ,  navgs,  f start,  fstop,  numpts) 


y,  Since  the  8722C  is  on  GPIB  address  7 
serial_write  (si ,  ’++addr  70 
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if  (nargin  >  2) 

7,  Ensure  proper  formatting  for  GPIB  command 
fstart  =  roundnCf start, -4) *10000; 
fstop  =  roundnCf stop, -4) *10000; 

numptsstr  =  [’POIN  ’  num2str(numpts) 
startstr  =  [’STAR  ’  num2str (fstart)  ’.E+5;’]; 
stopstr  -  [’STOP  ’  num2str (fstop)  ’.E+5;’]; 

serial_write(sl ,  numptsstr); 
serial_write(sl ,  startstr); 
serial_write(sl ,  stopstr); 

end 

if  (nargin  <  2) 

navgs  =  1 ;  end 

for  c  =  1 : navgs 

[F,M(:,c)]  =  get_8722c_sweep(sl) ; 

disp(  [’Trace  ’  num2str(c)  ’  out  of  ’  num2str (navgs)] ) ; 

end 

response  =  serial_write(sl ,  ’CONT;’); 
response  =  serial_write (si ,  ’++loc’); 


•/  »/  »/  »/  ®/  0/  0/  0/  ®/  ®/  ®/  V  V  V  V ' 
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o/fl  /o  /o  /o/o/o/o/o/o/o/o/o/o/9/o/«/a/o  /o 


8/  0/  0/  0/  0/  0/  0/  0/  0/  0/  0, 


7.  File; 

7o  Date; 

%  Author; 


serial_interactive .m 
28-Sep-2009 
amcc@mit . edu 


7. 

7o  Makes  a  giant  (~800k)  buffer  for  reading  serial  data  and  runs 
%  a  continuous  loop  interacting  with  the  serial  port 


maxPts  =  1601; 
bytesperline  =  512; 

warning  off  serial ; fscanf ;unsuccessfulRead 
if  "isempty (instrf ind) 

f close(instrfind) ;  end 
si  =  seriaK ’COM3' , ’BaudRate' ,115200,  ... 

’ InputBuf f erSize ’ ,maxPts*bytesperline , ’TimeOut ’,0.2); 
f open(sl) 

while  1 

s linput  =  input ( ’ Input  text ;  ’  ,  ’  s ’ ) ; 
response  =  serial_write(sl,  s 1 input ) ; 
dispC  [’response ;  ’  response]) 

end 


/e  /O  /9  /o  /O  /fl  /#  /O  /8  /O  /«  /9 


0/0/ 0/0/ 0/0/ 


'o  /o  /(I  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


0/  0/  0/  0/  0/  0/  0/  0/  0/  0/  0/ 


7»  File;  serial_readn.m 

7o  Date;  28-Sep-2009 
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%  Author:  amcc@mit.edu 

% 

%  Takes  in  a  serial  object  and  returns  n  lines  from  the  buffer 
y.  as  strings  in  a  list.  If  n  is  not  defined,  it  returns  one  line 
7,  not  in  a  list.  Terminates  if  the  line  received  is  empty 


function  response  =  serial_readn(sl ,  n) 

7.  Wait  for  the  transmission  to  finish 
BA  =  si .BytesAvailable; 
pause (0.5) 

while  (si. BytesAvailable  >  BA) 

BA  =  si .BytesAvailable ; 

disp( [’ [serial_readn]  Collecting  data...  Bytes  in  buffer: 
pause (1) 

end 

if  (nargin  <  2) 

response  =  fscanf  (si, ’y.s’)  ; 

else 

for  c  =  l:n 

response-[c}  =  fscanf  (si , ’7«s’ )  ; 
if  isempty(response{c}) ,  break,  end 

end 

end 


/o/o/o  /o/a/#/c  /•  /o 


/e  /o  /o  /o  /o  /9  /o  /e  /o  /•/o  /«  /o 


o/®/o/  0/  0/ 


num2str (BA)] ) ; 
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7.  File: 


serial_wait  .m 


7  Date:  28-Sep-2009 

7  Author:  ainccOmit.edu 

7 

7  Takes  in  a  serial  object  and  waits  up  to  n  seconds  for  a  single  message 
7  (ie  a  signal  message  that  some  process  is  done) 


fimction  response  =  serial_wait(sl ,  n) 

if  (nargin  <  2) 
n  =  1000;  end 

for  c  =  l:n 
pause ( . 2) 

response  =  fscanf (si, '7s’) ; 
if  "isempty (response) ,  break,  end 

disp( [’ [serial_wait]  No  response,  waiting  Is...  ’  num2str(c)  ’  elapsed’]) 
pause ( .8) 

end 


©/  #/  0/  0/  oy  0/  0/  oy  oy  oy  ft/  oy  ft/  ft/  oy  ft/  ft/  ft/  o/  o/  oy  «/  o/  «i 


/ft  /ft  /O  /ft  /ft  /ft 


0/  0/  0/  oy  oy  oy  oy  oy  «/  oy  oy  oy  oy  oy  oy  oy  oy  oy  o. 


^0  /ft  /o  /ft  /ft  /ft  /ft  /ft  /ft  /ft  /ft  /ft  /ft  /ft  /o  /ft  /ft  /ft  /ft  /•  /ft  /ft  /ft /ft  /ft  /ft  /ft  /«  /ft  /o 


7  File:  serial_write .m 

7  Date:  28-Sep-2009 

7  Author:  amcc@mit.edu 
7 

7  Writes  a  string  to  a  serial  object.  If  an  output  argximent  is 
7  defined,  returns  the  first  line  of  response 
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function  response  =  serial _write(sl ,  slinput) 
fprintf (si , slinput) 
if  nargout  >  0 
pause (0. 1) 

response  =  serial_readn(sl) ;  end 


end 


6/0/  0/  0/0/ 0/6/ 


/o/o/o  /e  ft 


6/  0/  0/  0/  6/  0/  0/  9/  6/  0/  6/  0/  0/  9/  0/  0/  0/  0/ 


°/o  File:  slist2nums .m 

7,  Date:  29-Sep-2009 

7o  Author :  amccSmit .  edu 


7. 

7.  Takes  in  a  list  of  strings,  reads  them  according  to  form  and  outputs 
7.  the  result 


function  nums  =  slist2nums(slist,  form) 

for  c  =  1 : length (slist) 
s  =  slist{c}; 

nums(c,:)  =  textscan(s,  form,  ’CollectOutput ’ , 1) ; 

end 

nums  =  cell2mat (nums)  ; 


/O  ft  ft  ft  ft  ft  ft  ft  ft 


0/  0/  0/  0/  0/  0/  0/  0/  0/  0/  0/  0/  9/  9/  0/  0/  9/  6/  6, 
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Appendix  B 


MATLAB  Surface  Impedance 
Numerical  Solver  Code 


0  /o  /o  h  /#  /o  /o  /o  /o  /o  /o  /o  /&  /«  /«  /o  /o  /o  /o 


function  [Zs  qq  K  Y]  =  surimp6(M,T,F,n\imqpts) 


y,  MAYBE  Write  out  exponential  form  for  Regime  III  to  speed  things  up 
y,  TODO:  Change  Zsint  to  eqn  B.8  pl48  Gao  to  remove  singularity 


uO  =  4*pi*le-7; 
w  =  2*pi*F; 

y.  Regime  I:  q  «  (max(  [1/M.xiO  1/M.l_mfp])) 

%  Should  be  constant  for  q  «  5e7 

y,  Regime  III:  q  »  (max(  [1/M.xiO  1/M.l_mfp])) 

y.  Should  be  constant  for  q  »  5e7 
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figure(82);  close(82);  figure(82); 

qregiiiie2  =  max([l/M.xiO 

KO  =  mbkernel6(qregime2/2e4,  w,  T,  M) ; 

YO  =  log(l+K0./(qregime2/(2e4)) ."2) ; 
semilogx(qregime2/2e4,Y0, ’bo’) ; 

global  qstore  Kstore  Ystore 

qstore  =  □  ;  Kstore  =  []  ;  Ystore  =  []  ; 

maxint  =  100;  rtol  =  le-6; 

QY  =  quadgk(®(q)mbkloop(q,w,T,M) ,  0,  Inf , ’MaxIntervalCount ’, maxint , ’RelTol ’, rtol , ’AbsTol ’ 

Zs  =  li*uO*w=t=pi/QY; 
qq  =  qstore; 

K  =  Kstore; 

Y  =  Ystore; 

function  Y  =  mbkloop(qq,w,T,M) 
global  qstore  Kstore  Ystore 


qregime2  =  max([l/M.xiO  1/M.l_mfp]); 

KO  =  mbkernel6(qregime2/2e4,  w,  T,  M) ; 

figure (82) 

K  =  zeros (size (qq) ) ; 

Y  =  zeros (size (qq) ) ; 
for  c  =  l:length(qq) 


if  qq(c)  <  qregime2/2e4 
K(c)  =  KO; 
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Y(c)  =  log(l+KO./qq(c)  .'■2); 


else 

K(c)  =  mbkerneie (qq(c) ,  w,  T,  M) ; 

Y(c)  =  log(l+K(c) ./qq(c) .~2)  ; 

end 

y,  disp([’q  =  ’  num2str(qq(c) )  K  =  ’  nuin2str(K(c) )] ) 
hold  on;  semilogx(qq(c) ,Y(c) , ’r . ’ ) ;  hold  off 

end 

qstore  =  [qstore  qq]  ; 

Kstore  =  [Kstore  K]  ; 

Ystore  =  [Ystore  Y] ; 

[qstore  qi]  =  sort (qstore) ; 

Kstore  =  Kstore (qi); 

Ystore  =  Ystore (qi); 


/o  /e  /o  /o  /o  /o  /o 


/o  /o/o/o/#  /o 


function  out  =  mbkernel6(q,  w,  T,  M) 
heV  =  4.135e-15;  7.  Planck  constant  [eV*s] 

hbareV  =  heV/(2*pi);  7.  Rationalized  planck  constant  [eV*s] 

[hw  D  D1  D2  b  al  a2  ap  am  gE  fE  fEhw]  =  mbkparams(l , 1 ,w,T,M) ; 

if  size(q,2)  >  size(q,l) 
q  =  transpose (q) ;  end 

rtol  =  le-5; 
maxint  =  100600; 
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y,  inf  in  =  D*5; 
7,  d  =  D=t‘rtol; 


7,  E  =  linspace(mcix(  [(D-hw+d)  -D]  )  ,  D-d,  numpts)  ; 

7.  Q1  =  trapz(E,Krealintl(E,q,w,T,M))  ; 

7.  [EM  qM]  =  meshgrid(E,q) ; 

7.  Q1  =  trapz(E,Krealintl(EM,qM,w,T,M)  ,2)  ; 

[Q1  errl]  =  quadgk(®(E)Krealintl (E,q,w,T,M) ,  max([(D-hw)  -D] ) ,  D, ’MaxIntervalCount’ .maxin 

if  (hw  <  2*D) 

Q2  =  zeros(size(Ql)) ; 

else 

7,  E  =  linspace(D-hw,  -D,  numpts); 

7,  [EM  qM]  =  meshgrid(E,q)  ; 

7,  Q2  =  trapz(EM,Krealint2(EM,qM,w,T,M)  ,2)  ; 

7.  Q2  =  trapz(E,Krealint2(E,q,w,T,M))  ; 

[Q2  err2]  =  quadgk(@(E)Krealint2(E,q,w,T,M) ,  D-hw,  -D, ’MaxIntervalCount \maxint , ’RelT 


end 

[Q3  err3]  =  quadgk(@(E)Krealint3(E,q,w,T,M) ,  D,  Inf , ’MaxIntervalCount’ ,maxint, ’RelTol’ ,rt 
[Q4  err4]  =  quadgk(@(E)Krealint4(E,q,w,T,M) ,  D,  Inf , 'MaxIntervalCount ’ .maxint , ’RelTol’ ,rt 

7,  E  =  logspace(loglO(D+d)  ,  loglO(D+d)  +  loglOCl  .0001)  ,  numpts); 

7.  E  =  [E  logspace(loglO(D+d)  +  loglOCl  .0001)  ,  logl0(inf  in) ,  numpts)]; 

7,  [EM  qM]  =  meshgrid(E,q)  ; 

7,  Q3  =  trapz(E,Krealint3(EM,qM, w,T,M)  , 2)  ; 

7.  Q4  =  trapz(E,Krealint4(EM,qM,w,T,M)  ,2)  ; 

7«  Q3  =  trapz(E,Krealint3(E,q,w,T,M))  ; 

7.  Q4  =  trapz(E,Krealint4(E,q,w,T,M))  ; 
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if  (hw  <  2=t=D) 

Q5  =  zeros (size (Ql) ) ; 

else 

7.  E  =  linspace(D-hw,  -D,  numpts)  ; 

7.  [EM  qM]  =  nieshgrid(E,q)  ; 

7,  Q5  =  trapz(E,Kimagint5(EM,qM,w,T,M)  ,2) ; 

7.  Q5  =  trapz(E,Kimagint5(E,q,w,T,M))  ; 

[q5  err5]  =  quadgk(@(E)Kimagint5(E,q,w,T,M),  max([(D-hw)  -D] ) ,  -D, ’MaxIntervalCount ’ 

end 

%  E  =  linspace(D+d,  infin,  numpts); 

7.  [EM  qM]  =  meshgrid(E,q) ; 

7»  Q6  =  trapz(E,Kimagint6(EM,qM,w,T,M)  ,2)  ; 

%  Q6  =  trapz(E,Kimagint6(E,q,w,T,M)) ; 

[Q6  err6]  =  quadgk(®(E)Kimagint6(E,q,w,T,M) ,  D,  Inf , ’MaxIntervalCount'.maxint, 'RelTol’ ,rt 


7. 

7. 

7,  dispC  [’Integral  1)  ’  num2str(Ql)  err  of  ’  num2str(errl/Ql*100)  ’7o’]) 
7o  disp(  [’Integral  2)  ’  num2str(Q2)  ’,  err  of  ’  num2str (errl/Q2*100)  ’7.’]) 

7,  disp(  [’Integral  3)  ’  num2str(Q3)  ’,  err  of  ’  num2str(errl/Q3*100)  ’7.’]) 

7.  disp([’ Integral  4)  ’  num2str(Q4)  ’,  err  of  ’  num2str (errl/Q4*100)  ’7.’]) 
7,  disp([’ Integral  5)  ’  num2str(Q5)  ’,  err  of  ’  num2str (errl/QS^lOO)  ’7.’]) 

7o  disp(  [’Integral  6)  ’  num2str(Q6)  ’,  err  of  ’  num2str(errl/Q6*100)  ’7.’]) 


Kreal  =  3./(q*libareV*M.nu0*M.L0''2)  .*(Q1+1  ./2.*Q2-Q3+Q4)  ; 
Kimag  =  3./(q*bbareV*M.nu0*M.L0~2) .*(-l./2.*Q5+Q6) ; 


out  =  transpose (Kreal  +  li. ♦Kimag) ; 
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function  plotdebug(E,Y,n) 

figure(n) ;  plot(E,abs(Y) , ’b-’ ,E,real(Y) , ’x-’ ,E,imag(Y) , ’g-’) ; 

7,  disp([’Sum  ’  num2str(n)  ’)  real:  ’  num2str (sum(real (Y) ) )  imag:  '  num2str (sum(imag(Y) 

%  global  EE  YY 
7.  EE  =  E; 

%  YY  =  Y; 

%  pause 

function  Y  =  Krealintl(E,q,w,T,M) 

[hw  D  D1  D2  b  al  a2  ap  am  gE  fE  fEhw]  =  mbkparams(E,q,w,T,M) ; 

Y  =  (l-2.*fEhw)  . . . 

((E.~2+D~2+hw*E) ./(sqrt(D~2-E.‘2) .^squbCCE+hw) ."2-D''2))  .*  R(a2,al+b)  +  S(a2,al+b) 
7.  Y  =  real(Y);  7.  Since  it’s  an  integral  for  the  real  part 
7o  plotdebug(E,Y,  1) ; 

fimction  Y  =  Krealint2(E,q,w,T,M) 

[hw  D  D1  D2  b  al  a2  ap  am  gE  fE  fEhw]  =  mbkparams(E,q,w,T,M) ; 

Y=  (l-2.*fEhw) .*((gE+l) .*S(am.b)  -  (gE-1) .*S(ap,b)) ; 

7»  plotdebug(E,Y,2)  ; 


function  Y  =  Krealint3(E,q,w,T,M) 

[hw  D  D1  D2  b  al  a2  ap  am  gE  fE  fEhw]  =  mbkparams(E,q,w,T,M) ; 

Y  =  (1-fE-fEhw) .*(gE-l) .*S(ap,b) ; 

7o  plotdebug(E,Y,3)  ; 
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function  Y  =  Krealint4(E,q,w,T,M) 

[hw  D  D1  D2  b  al  a2  ap  am  gE  fE  fEhw]  =  mbkparams(E,q,w,T,M) ; 

Y  =  (fE  -  fEhw) .*(gE+l) .*S(am,b) ; 

“/o  plotdebug(E,Y,4)  ; 


fiinction  Y  -  Kimagint5(E,q,w,T,M) 

[hw  D  D1  D2  b  al  a2  ap  am  gE  fE  fEhw]  =  mbkparams(E,q,w,T,M) ; 

Y  =  (l-2.*fEhw) .*((gE+l) .*R(am,b)+(gE-l) .*R(ap,b)) ; 

%  plotdebugCE, Y, 5) ; 


function  Y  =  Kimagint6(E,q,w,T,M) 

[hw  D  D1  D2  b  al  a2  ap  am  gE  fE  fEhw]  =  mbkparams(E,q,w,T,M) ; 

Y  =  (fE-fEhw) .*((gE+l) .*R(am,b)+(gE-l) .*R(ap,b)) ; 

%  plotdebug(E,Y,6) ; 


function  out  =  R(a,b) 
7,  Updated  5/19/2010 


s  =  b-li*a; 

out  =  real(-s/2+(s. ~2+l)/2.*atan(l./s)) ; 

7o  nx  =  zeros  (size  (a) ) ; 

7,  nx(b.  "2+a.  "2-1  <  0)  =  ones(size(nx(b."2+a."2-l  <  0))); 

7. 
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7o  out  =  -b./2  +  a. +b . /4 .  ^logC  (b .  ~2+(l+a)  .  ”2)  . /(b . ''2+(l-a) .  ”2))  ... 

%  +  l./4.*(l+b.''2-a.~2)  .=t'(atan(2.*b./(b."2+a.''2-l))  +nx.*pi); 


%  Eqn  A. 25  from  Appendix  A,  Gao  thesis  pl45 
function  out  =  S(a,b) 

7,  Updated  5/19/2010 


s  =  b-li*a; 

out  =  imag(-s/2+(s  . ''2+l)/2.  *atan(l./s) )  ; 

7.  nx  =  zeros  (size  (a)  )  ; 

7.  nx(b. ''2+a.  “2-1  <  0)  =  ones(size(nx(b.  “2+a.  “2-1  <  0))); 

7. 

7.  out  =  a./2  -  a.=i‘b./2.=i‘(atan(2.*b./(b.“2+a.“2-l))+nx.*pi)  ... 

7o  +  1 .  /8.*(l+b.  “2-a.  “2)  .*log((b.  “2+(l+a)  .  “2)  ./(b.  “2+(l-a)  .  “2))  ; 

function  fE  =  fermi(E,T) 
kbeV  =  8.617e-5;  7.  eV./K 


if  (T==0) 

fE  =  zeros (size (E) ) ; 

fE(E<=0)  =  ones(size(E(E<=0))) ; 

else 

fE  =  l./(l+exp(E/(kbeV*T))) ; 

end 


function  gE  =  g(E,hw,D,Dl,D2) 


89 


gE  =  (E.''2+D.''2+hw.*E)  ./(D1.+D2)  ; 


fxinction  [hw  D  D1  D2  b  al  a2  ap  am  gE  fE  fEbw]  =  mbkparams(E,q,w,T,M) 

heV  =  4.135e-15;  7.  Planck  constant  [eV*s] 

hbareV  =  heV/(2*pi);  %  Rationalized  planck  constant  [eV*s] 

hw  =  hbareV  *  w; 
b  =  1 . / (q*M. l_mfp) ; 

D  =  reducedenergygap(T,M.Tc,M.DO) ; 

D1  =  sqrt (E. ~2-D . ~2) ; 

D2  =  sqrt ( (E+hw) . ~2-D. “2) ; 
al  =  D1 . / (hbareV*M .nuO*q) ; 
a2  =  D2./(hbareV*M.nuO*q) ; 
ap  =  al+a2; 
am  =  a2-al; 
gE  =  g(E,hw,D,Dl,D2) ; 
fE  =  fermi (E,T) ; 
fEhw  =  fermi (E+hw,T) ; 


y "/  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 

/«  /o/o  /o  /q  /«  /o  /o  /o  /o  /o  /e  /e  /q  /o  /o  /o  /«  /e  /o  /•  /o  /«  /o  /«  /o  /o  /o  In  h  In  in  in  In  in  in  In  in  in  in  in  in  in  in  in  in  in  in  in  in  in  in 


yyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyyy 

/o  It  ft  Itltltltltltltltltltltltltltlt  It  Itititititit  It  It  hit  It  It  h  h  Ititit  h  h  Ititit  h  h  h  Ititit  h  h  It  It  It  h 


7,  As  described  on  p24  Gao  thesis,  section  2.2.6 
7o  Table  data  taken  from  Muhlschlegel  for  T/Tc  >  0.18 
7.  Adapted  verbatim  from  SuperMix  code  in  supercond.cc 


function  D  =  reducedenergygap(T,Tc ,D0) 
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X  =  T/Tc; 
xmatch  =  0.32; 

ymatch  =  exp(-sqrt(3.562*xmatch)*exp(-1.764/xmatch)) ; 

y,  //  Table  of  reduced  energy  gap  vs.  T/Tc,  from  Muhlschlegel  (1959). 
y,  //  Table  starts  at  x  =  T/Tc  =  0.18,  and  goes  in  steps  of  dx  =  0.02 
ratiotable  =  [1.0,  0.9999,  0.9997,  0.9994,  0.9989,  ... 

0.9982,  0.9971,  0.9957,  0.9938,  0.9915,  ... 

0.9885,  0.985,  0.9809,  0.976,  0.9704,  ... 

0.9641,  0.9569,  0.9488,  0.9399,  0.9299,  ... 

0.919,  0.907,  0.8939,  0.8796,  0.864,  ... 

0.8471,  0.8288,  0.8089,  0.7874,  0.764,  ... 

0.7386,  0.711,  0.681,  0.648,  0.6117,  ... 

0.5715,  0.5263,  0.4749,  0.4148,  0.3416,  ... 

0.2436,  0.0]; 

y,  static  double  gap_supcond (double  x) 

%  -c 

y,  //  the  argument  x  =  T/Tc;  check  if  x  is  physical 

y,  if(x  <  0.  II  X  >=  1.)  return(O.O);  //  return  gap  of  zero 

if  (x>l)  II  (x<0) ,  error ( ’Nonphysical  x’);  end 

y. 

y,  //  if  T/Tc  <  0.04,  reduced  energy  gap  is  1.0,  to  better  than  1  part  in  10“22 

y,  if(x  <  0.04)  return (1.0);  //  return  gap  of  1.0 

if  (x<0.04),  ratio  =  1; 

y,  //  if  T/Tc  <=  xmatch,  we  use  an  analytic  expression  to  extend  the  gap  table; 
"/o  //  (  which  uses  k  Tc  “=  Delta(0)/1.764;  c.f.  Tinkham  ) 
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7. 


7,  static  const  double  xmatch  =  0.32;  //  xmatch  must  be  a  multiple  of  0.02,  >=  0.18 

7o  static  const  double  ymatch  =  exp(-sqrt(3.562*xmatch)*exp(-1.764/xmatch)) ; 

7o  if(x  <=  xmatch)  return(exp(-sqrt(3.562+x)*exp(-1.764/x))) ;  // 

7.  3.562==2*Pi/1.764 

elseif  (x  <  xmatch) 

ratio  =  exp(-sqrt(2*pi/1.764=t=x)=t=exp(-1.764/x))  ; 

%  //  if  we  get  here,  we  interpolate  between  data  points  to  calculate  gap,  taking 

7o  //  sqrtd  -  T/Tc)  dependence  near  the  gap  into  account. 

7. 

%  int  index  =  static_cast<int>(floor((x-0. 18)/0.02)) ;  //  index  into  ratio [] 

7o  double  xl  =  index*0.02  +  0.18, 

7»  XU  =  xl  +  0.02, 

7o  yl  =  ((xl  >  xmatch)  ?  ratio[index]  ;  ymatch)/sqrt(l.-xl) , 

7o  yu  =  (index  <  40)  ?  ratio  [index+1] /sqrt(l . -xu)  :  1.74; 

7o  //  Behavior  near  Tc  -  see  Tinkham,  eqn.  2-54. 

else 

index  =  floor((x-0. 18)/0.02) ; 

xl  =  (index-1) *0.02  +  0.18;  7o  index-1  since  in  C++  first  point  is  0 
XU  =  xl  +  0.02; 

if  (xl  >  xmatch) 

yl  =  ratiotable (index) /sqrt ( 1-xl) ; 

else 

yl  =  ymatch/sqrt(l-xl) ; 

end 

if  (index  <  40) 

yu  =  ratiotable (index+1) /sqrt ( 1-xu) ; 


92 


else 


yu  =  1.74; 

end 


ratio  =  sqrtCl  -  x)  *  (yl  +  (yu-yD  +  Cx-xD/Cxu-xl)) 

end 


D  =  DO*ratio; 


's/o/o/o/o/oA/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/oA/o/e/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/o/fl 
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